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Preface 


The  study  documented  herein  was  undertaken  as  part  of  Work  Unit  387- 
9456h,  “Seismic  Design  of  Cantilever  Retaining  Walls,”  funded  by  the  Head¬ 
quarters,  U.S.  Army  Corps  of  Engineers  (HQUSACE)  Civil  Works  Earthquake 
Engineering  Research  Program  (EQEN)  under  the  purview  of  the  Geotechnical 
and  Structures  Laboratory  (GSL),  Vicksburg,  MS,  U.S.  Army  Engineer  Research 
and  Development  Center  (ERDC).  Technical  Director  for  this  research  area  was 
Dr.  Mary  Ellen  Hynes,  GSL.  The  HQUSACE  Program  Monitor  for  this  work 
was  Ms.  Anjana  Chudgar.  The  principal  investigator  (PI)  for  this  study  was 
Dr.  Robert  M.  Ebeling,  Computer-Aided  Engineering  Division  (CAED),  Infor¬ 
mation  Technology  Laboratory  (ITL),  Vicksburg,  MS,  ERDC,  and  Program 
Manager  was  Mr.  Donald  E.  Yule,  GSL.  The  work  was  performed  at  University 
of  Michigan,  Ann  Arbor,  and  at  ITL.  The  effort  at  the  University  of  Michigan 
was  funded  through  response  to  the  ERDC  Broad  Agency  Announcement  FYOl, 
BAA#  ITL- 1,  “A  Research  Investigation  of  Dynamic  Earth  Loads  on  Cantilever 
Retaining  Walts  as  a  Function  of  the  Wall  Geometry,  Backfill  Characteristics, 
and  Numerical  Modeling  Technique.” 

This  research  was  performed  and  the  report  prepared  by  Dr.  Russell  A.  Green 
of  the  Department  of  Civil  and  Environmental  Engineering,  University  of 
Michigan,  and  by  Dr.  Ebeling  under  the  direct  supervision  of  Mr.  H.  Wayne 
Jones,  CAED,  and  Dr.  Jeffery  P.  Holland,  Director,  ITL.  The  work  was 
performed  during  the  period  December  2001  to  August  2002  by  Dr.  Green  and 
Dr.  Ebeling.  This  report  summarizes  the  results  of  the  first  phase  of  a  research 
investigation  examining  the  seismic  loads  induced  on  the  stem  of  a  cantilever 
retaining  wall.  This  investigation  marks  the  first  use  of  the  computer  program 
FLAG  (Fast  Lagrangian  Analysis  of  Continua)  for  analyzing  the  dynamic 
response  of  a  Corps  earth  retaining  structure,  with  the  emphasis  of  the 
investigation  being  on  the  details  of  numerical  modeling  with  FLAC,  as  well  as 
the  results  of  the  analyses.  Further  analyses  are  required  to  confirm  the  identified 
trends  in  the  results  of  the  analyses  and  to  formulate  design  recommendations  for 
Corps  earth  retaining  structures.  During  the  course  of  this  research  investigation, 
the  authors  had  numerous  discussions  with  other  FLAC  users.  Of  particular  note 
were  the  lengthy  conversations  with  Mr.  Guney  Olgun,  Virginia  Polytechnic  and 
State  University,  Blacksburg,  which  were  instrumental  in  completing  Phase  1  of 
this  research  investigation.  Others  who  provided  valuable  insight  into  the 
workings  of  FLAC  were  Mr.  Nason  McCullough  and  Dr.  Stephen  Dickenson, 
Oregon  State  University,  Corvallis;  Dr.  N.  Deng  and  Dr.  Farhang  Ostadan, 
Bechtel  Corporation,  San  Francisco,  CA;  Mr.  Michael  R.  Lewis,  Bechtel 


Savannah  River,  Inc.,  Aiken,  SC;  Dr.  Peter  Byrne  and  Dr.  Mike  Beaty, 
University  of  British  Columbia,  Vancouver;  and  Dr.  Marte  Gutierrez,  Virginia 
Tech. 

At  the  time  of  publication  of  this  report.  Dr.  James  R.  Houston  was  Director, 
ERDC,  and  COL  John  W.  Morris  III,  EN,  was  Commander  and  Executive 
Director. 


The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication, 
or  promotional  purposes.  Citation  of  trade  names  does  not  constitute  an 
official  endorsement  or  approval  of  the  use  of  such  commercial  products. 
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Introduction 


1.1  Introduction 

This  report  presents  the  results  of  the  first  phase  of  a  research  investigation 
into  the  seismic  response  of  earth  retaining  structures  and  the  extension  of  the 
displacement  controlled  design  procedure,  as  applied  to  the  global  stability 
assessment  of  Corps  retaining  structures,  to  issues  pertaining  to  their  internal 
stability.  It  is  intended  to  provide  detailed  information  leading  to  refinement  of 
the  Ebeling  and  Morrison  (1992)  simplified  seismic  engineering  procedure  for 
Corps  retaining  structures.  Specific  items  addressed  in  this  Phase  1  report  deal 
with  the  seismic  loads  acting  on  the  stem  portion  of  cantilever  retaining  walls.  A 
typical  Corps  cantilever  retaining  wall  is  shown  in  Figure  1-1.  It  is  envisioned 
that  this  information  will  be  used  in  the  development  of  a  refined  engineering 
procedure  of  the  stem  and  base  reinforced  concrete  cantilever  wall  structural 
members  for  seismic  structural  design. 


Figure  1-1.  Typical  Corps  cantilever  wall,  including  structural  and  driving  wedges 
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1.2  Background 

Formal  consideration  of  the  permanent  seismic  wall  displacement  in  the 
seismic  design  process  for  Corps-type  retaining  structures  is  given  in  Ebeling  and 
Morrison  (1992).  The  key  aspect  of  this  engineering  approach  is  that  simplified 
procedures  for  computing  the  seismically  induced  earth  loads  on  retaining 
structures  are  dependent  upon  the  amount  of  permanent  wall  displacement  that  is 
expected  to  occur  for  each  specified  design  earthquake.  The  Corps  uses  two 
design  earthquakes  as  stipulated  in  Engineer  Regulation  (ER)  1 1 10-2-1806 
(Headquarters,  U.S.  Army  Corps  of  Engineers  (HQUSACE)  1995):  the 
Operational  Basis  Earthquake  (OBE)'  and  the  Maximum  Design  Earthquake 
(MDE).  The  retaining  wall  would  be  analyzed  for  each  design  case.  The  load 
factors  used  in  the  design  of  reinforced  concrete  hydraulic  structures  are  different 
for  each  of  these  two  load  cases. 

The  Ebeling  and  Morrison  simplified  engineering  procedures  for  Corps 
retaining  structures,  as  described  in  their  1992  report,  are  geared  toward  hand 
calculations.  However,  research  efforts  are  currently  underway  at  the  U.S.  Army 
Engineer  Research  and  Development  Center  (ERDC)  to  computerize  these 
engineering  procedures  and  to  make  possible  the  use  of  acceleration  time- 
histories  in  these  design/analysis  processes  when  time-histories  are  made 
available  on  Corps  projects.  In  the  Ebeling  and  Morrison  simplified  seismic 
analysis  procedure  two  limit  states  are  established  for  the  backfill;  the  first 
corresponds  to  walls  retaining  yielding  backfill,  while  the  second  corresponds  to 
walls  retaining  nonyielding  backfill.  Examples  of  Corps  retaining  walls  that 
typically  exhibit  these  two  conditions  in  seismic  evaluations  are  shown  in  Fig¬ 
ure  1-2.  In  this  figure  Fv  and  F  h  are  the  vertical  and  horizontal  components, 
respectively,  of  the  resultant  force  of  the  stresses  acting  on  imaginary  sections 
A-A  and  B-B,  and  Tand  N'  are  the  shear  and  normal  reaction  forces,  respectively, 
on  the  bases  of  the  walls. 

It  is  not  uncommon  for  retaining  walls  of  the  type  shown  in  Figure  l-2a,  i.e., 
soil-founded  cantilever  retaining  walls,  to  have  sufficient  wall  movement  away 
from  the  backfill  during  a  seismic  event  to  mobilize  the  shear  strength  within  the 
backfill,  resulting  in  active  earth  pressures  acting  on  the  structural  wedge  (as 
delineated  from  the  driving  wedge  by  imaginary  section  A-A  extending  vertically 
from  the  heel  of  the  wall  up  through  the  backfill).  Figure  l-2b  shows  a  wall 
exemplifying  the  second  category,  walls  retaining  a  nonyielding  backfill.  For  a 
massive  concrete  gravity  lock  wall  founded  on  competent  rock  with  high  base 
interface  and  rock  foundation  shear  strengths  (including  high-  strength  rock 
joints,  if  present,  within  the  foundation),  it  is  not  uncommon  to  find  that  the 
typical  response  of  the  wall  during  seismic  shaking  is  the  lock  wall  rocking  upon 
its  base.  For  this  case,  wall  movements  in  sliding  are  typically  not  sufficient  to 
mobilize  the  shear  strength  in  the  backfill. 
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‘  For  convenience,  symbols  and  unusual  abbreviations  are  listed  and  defined  in  the 
Notation  (Appendix  E). 
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Figure  1-2.  Earth  retaining  structures  typical  of  Corps  projects;  (a)  soil-founded, 
cantilever  floodwall  retaining  earthen  backfill;  (b)  rock-founded, 
massive  concrete  lock  wall  retaining  earthen  backfill 

Yielding  backfills  assume  that  the  shear  strength  of  the  backfill  is  fully 
mobilized  (as  a  result  of  the  wall  moving  away  from  the  backfill  during  earth¬ 
quake  shaking),  and  the  use  of  seismically  induced  active  earth  pressure  relation¬ 
ships  (e.g.,  Mononobe-Okabe)  is  appropriate.  A  calculation  procedure  first 
proposed  by  Richards  and  Elms  (1979)  for  walls  retaining  “dry”  backfills  (i.e.,  no 
water  table)  is  used  for  this  limit  state.  Ebeling  and  Morrison  (1992)  proposed 
engineering  calculation  procedures  for  “wet”  sites  (i.e.,  sites  with  partially  sub¬ 
merged  backfills  and  for  pools  of  standing  water  in  the  chamber  or  channel)  and 
developed  a  procedure  to  compute  the  resultant  active  earth  pressure  force  acting 
on  the  structural  wedge  using  the  Mononobe-Okabe  relationship.  (Most  Corps 
sites  are  “wet”  since  the  Corps  usually  deals  with  hydraulic  structures.)  The 
simplified  Ebeling  and  Morrison  engineering  procedure  recommends  that  a 
Richards  and  Elms  type  displacement-controlled  approach  be  applied  to  the  earth 
retaining  structure,  as  described  in  Section  6.3  of  Ebeling  and  Morrison  (1992) 
for  Corps  retaining  structures.  It  is  critical  to  the  calculations  that  partial  sub¬ 
mergence  of  the  backfill  and  a  standing  pool  of  water  in  the  chamber  (or  channel) 
are  explicitly  considered  in  the  analysis,  as  given  by  the  Ebeling  and  Morrison 
simplified  computational  procedure.  Equations  developed  by  Ebeling  and  Morri¬ 
son  to  account  for  partial  submergence  of  the  backfill  in  the  Mononobe-Okabe 
resultant  active  earth  pressure  force  computation  is  given  in  Chapter  4  of  their 
report.  A  procedure  for  assigning  the  corresponding  earth  pressure  distribution 
was  developed  by  Ebeling  and  Morrison  for  a  partially  submerged  backfill  and  is 
described  using  Figures  7.8,  7.9,  and  7.10  of  their  report. 

Key  to  the  categorization  of  walls  retaining  yielding  backfills  in  the  Ebeling 
and  Morrison  simplified  engineering  procedure  for  Corps  retaining  structures  is 
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the  assessment  by  the  design  engineer  of  the  minimum  seismically  induced  wall 
displacements  to  allow  for  the  full  mobilization  of  the  shear  resistance  of  the 
backfill  and,  thus,  the  appropriate  use  of  the  Mononobe-Okabe  active  earth 
pressure  relationship  in  the  computations.  Ebeling  and  Morrison  made  a  careful 
assessment  of  the  instrumented  dynamic  earth  pressure  experiments  available  in 
the  technical  literature  prior  to  their  publication  in  1992.  The  results  of  this 
assessment  are  described  in  Chapter  2  of  Ebeling  and  Morrison  (1992).  Ebeling 
and  Morrison  concluded  that  the  minimum  wall  displacement  criteria  developed 
by  Clough  and  Duncan  (1991)  for  the  development  of  “active”  static  earth 
pressure  are  also  reasonable  guidance  for  the  development  of  seismically  induced 
active  earth  pressure.  This  guidance  for  engineered  backfills  is  given  in  Table  1 
of  Ebeling  and  Morrison  (1992).  Minimum  permanent  seismically  induced  wall 
displacements  away  from  the  backfill  are  expressed  in  this  table  as  a  fraction  of 
the  height  of  backfill  being  retained  by  the  wall.  The  value  for  this  ratio  is  also  a 
function  of  the  relative  density  of  the  engineered  backfill.  Thus,  prior  to 
accepting  a  permanent  seismic  wall  displacement  prediction  made  following  the 
simplified  displacement-controlled  approach  for  Corps  retaining  structures 
(Section  6.3  of  Ebeling  and  Morrison  1992),  the  design  engineer  is  to  check  if  his 
computed  permanent  seismic  wall  displacement  value  meets  or  exceeds  the 
minimum  displacement  value  for  active  earth  pressure  given  in  Table  1  of 
Ebeling  and  Morrison  (1992).  This  ensures  that  the  use  of  active  earth  pressures 
in  the  computation  procedure  is  appropriate. 

In  the  second  category  of  walls  retaining  nonyielding  backfills  (Figure  l-2b), 
Ebeling  and  Morrison  recommend  the  use  of  at-rest  type,  earth  pressure 
relationship  in  the  simplified  hand  calculations.  Wood's  (1973)  procedure  is  used 
to  compute  the  incremental  pseudo-static  seismic  loading,  which  is  superimposed 
on  the  static,  at-rest  distribution  of  earth  pressures.  Wood’s  is  an  expedient  but 
conservative  computational  procedure  (Ebeling  and  Morrison  (1992),  Chapter  5). 
(A  procedure  to  account  for  wet  sites  with  partially  submerged  backfills  and  for 
pools  of  standing  water  in  the  chamber  or  channel  was  developed  by  Ebeling  and 
Morrison  (1992)  and  outlined  in  Chapter  8  of  their  report.)  It  is  Ebeling’s 
experience  with  the  type  lock  walls  shown  in  Figure  l-2b  of  dimensions  that  are 
typical  for  Corps  locks  that  seismically  induced  sliding  is  an  issue  only  with  large 
ground  motion  design  events  and/or  when  a  weak  rock  joint  or  a  poor  lock-to- 
foundation  interface  is  present. 

After  careful  deliberation,  Ebeling  and  Morrison  in  consultation  with  Whit- 
man‘  and  Finn^  judged  the  simplified  engineering  procedure  for  walls  retaining 
nonyielding  backfills  applicable  to  walls  in  which  the  wall  movements  are  small, 
less  than  one-fourth  to  one-half  of  the  Table  1  (Ebeling  and  Morrison  1992) 
active  displacement  values.  Recall  that  the  Ebeling  and  Morrison  engineering 
procedure  is  centered  on  the  use  of  one  of  only  two  simplified  hand- 
computational  procedures. 


*  Dr.  Robert  V.  Whitman,  1992,  Professor  Emeritus,  Massachusetts  Institute  of 
Technology,  Boston. 

^  Dr.  W.  D.  Liam  Finn,  1992,  Professor  Emeritus,  University  of  British  Columbia, 
Vancouver. 
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Rotational  response  of  the  wall  (compared  to  sliding)  is  beyond  the  scope  of 
the  Ebeling  and  Morrison  (1992)  simplified  engineering  procedures  for  Corps 
retaining  structures.  This  1992  pioneering  effort  for  the  Corps  dealt  only  with  the 
sliding  mode  of  permanent  displacement  during  seismic  design  events.  It  is 
recognized  that  the  Corps  has  some  retaining  structures  that  are  more  susceptible 
to  rotation-induced  (permanent)  displacement  during  seismic  events  than  to 
(permanent)  sliding  displacement.  To  address  this  issue,  Ebeling  is  currently 
conducting  research  at  ERDC  leading  to  the  development  of  a  simplified  engi¬ 
neering  design  procedure  for  the  analysis  of  retaining  structures  that  are  con¬ 
strained  to  rotate  about  the  toe  of  the  wall  during  seismic  design  events  (Ebeling 
emd  White,  in  preparation). 


1 .3  Research  Objective 

The  Ebeling  and  Morrison  (1992)  simplified  seismic  engineering  procedures 
for  Corps  retaining  structures  did  not  address  issues  pertaining  to  the  structural 
design  of  cantilever  retaining  walls.  The  objective  of  the  research  described  in 
this  report  is  to  fill  this  knowledge  gap  and  determine  the  magnitude  and  distribu¬ 
tion  of  the  seismic  loads  acting  on  cantilever  retaining  walls  for  use  in  the  design 
of  the  stem  and  base  reinforced  concrete  cantilever  wall  structural  members. 


1 .4  Research  into  the  Seismic  Response  of  a 
Cantilever  Retaining  Wall 

The  seismic  loads  acting  on  the  structural  wedge  of  a  cantilever  retaining 
wall  are  illustrated  in  Figure  1-3.  The  structural  wedge  consists  of  the  concrete 
wall  and  the  backfill  above  the  base  of  the  wall  (i.e.,  the  backfill  to  the  left  of  a 
vertical  section  through  the  heel  of  the  cantilever  wall).  The  resultant  force  of 
the  static  and  dynamic  stresses  acting  on  the  vertical  section  through  the  heel 
(i.e.,  heel  section)  is  designated  as  Pae,  heeh  and  the  normal  and  shear  base 
reactions  are  N'  and  T,  respectively.  Seismically  induced  active  earth  pressures 
on  the  heel  section,  Pae,  heeh  are  used  to  evaluate  the  global  stability  of  the 
structural  wedge  of  a  cantilever  retaining  wall,  presuming  there  is  sufficient  wall 
movement  away  from  the  backfill  to  fully  mobilize  the  shear  resistance  of  the 
retained  soil.  The  relative  slenderness  of  the  stem  portion  of  a  cantilever  wall 
requires  structural  design  consideration.  In  Figure  1-3  the  seismically  induced 
shear  and  bending  moments  on  a  section  of  the  stem  are  designated  as  s  and  w, 
respectively.  The  resultant  force  of  the  static  and  dynamic  stresses  acting  on  the 
stem  of  the  wall  shown  in  Figure  1-3  is  designated  as  Pe,  stem-  The  A  is  not 
included  in  the  subscript  because  the  structural  design  load  is  not  necessarily 
associated  with  active  earth  pressures. 

A  dry  site  (i.e.,  no  water  table)  will  be  analyzed  in  this  first  of  a  series  of 
analyses  of  cantilever  retaining  walls  using  FLAG  (Fast  Lagrangian  Analysis  of 
Continua).  This  allows  the  researchers  to  gain  a  full  understanding  of  the 
dynamic  behavior  of  the  simpler  case  of  a  cantilever  wall  retaining  dry  backfill 
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before  adding  the  additional  complexities  associated  with  submerged  or  partially 
submerged  backfills. 

This  report  summarizes  the  results  of  detailed  numerical  analyses  performed 
on  a  cantilever  wall  proportioned  and  structurally  detailed  per  Corps  guidelines 
given  in  Engineer  Manuals  (EM)  1 1 10-2-2104  (HQUSACE  1992)  and  1 1 10-2- 
2502  (HQUSACE  1989))  for  global  stability  and  structural  strength  under  static 
loading.  The  objective  of  the  analyses  was  to  identify  trends  and  correlations 
between  Pae.  hed  and  Pe.  sum  and  their  respective  points  of  application.  The  identi¬ 
fication  of  such  trends  allows  the  displacement-controlled  design  procedure, 
which  can  be  used  to  estimate  Pae.  hed,  to  be  extended  to  estimate  Pe.  stem,  which  is 
required  for  the  structural  design  of  the  stem. 

The  detailed  numerical  analyses  were  performed  using  the  commercially 
available  computer  program  FLAC.  The  nonlinear  constitutive  models,  in 
conjunction  with  the  explicit  solution  scheme,  in  FLAC  give  stable  solutions  to 
unstable  physical  processes,  such  as  the  sliding  or  overturning  of  a  retaining  wall. 
FLAC  allows  permanent  displacements  to  be  modeled,  which  is  inherently 
required  by  the  displacement-controlled  design  procedure.  The  resultant  forces 
acting  on  the  heel  sections  and  their  points  of  applications  as  determined  from  the 
FLAC  analyses  were  compared  with  values  computed  using  the  Mononobe- 
Okabe  equations  in  conjunction  with  the  displacement-controlled  design 
procedure  (e.g.,  Ebeling  and  Morrison  1992). 
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1.5  Organization  of  Report 

The  organization  of  the  report  follows  the  sequence  in  which  the  work  was 
performed.  Chapter  2  outlines  the  process  of  selecting  the  ground  motions  (e.g., 
acceleration  time-histories)  used  in  the  FLAG  analyses.  Chapter  3  gives  a  brief 
overview  of  the  numerical  algorithms  in  FLAG  and  outlines  how  the  various 
numerical  model  parameters  were  determined.  Chapter  4  describes  the  data 
reduction  and  interpretation  of  the  FLAG  results,  followed  by  the  References. 
Appendix  A  provides  detailed  calculation  of  the  geometry  and  structural  design 
for  static  loading  of  the  wall  analyzed  dynamically.  Appendix  B  reviews  the  sign 
convention  and  notation  used  in  this  report  and  also  presents  the  Mononobe- 
Okabe  earth  pressure  equations  (e.g.,  Ebeling  and  Morrison  1992,  Chapter  4). 
Appendix  C  is  a  brief  overview  of  the  displacement-controlled  procedure  for 
global  stability  of  retaining  walls.  Finally,  Appendix  D  summarizes  a  parameter 
study  performed  to  determine  how  best  to  specify  ground  motions  in  FLAG. 


1.6  Future  Work 

This  report  presents  the  results  of  the  first  phase  of  an  ongoing  research 
investigation.  Additional  FLAG  analyses  are  planned  to  determine  if  the 
observed  trends  presented  in  Chapter  4  of  this  report  are  limited  to  the  wall 
geometry  and  soil  conditions  analyzed,  or  whether  they  are  general  trends  that 
are  applicable  to  other  wall  geometries  and  soil  conditions.  Additionally,  the 
same  walls  analyzed  using  FLAG  will  be  analyzed  using  the  computer  program 
FLUSH.  FLUSH  solves  the  equations  of  motions  in  the  frequency  domain  and 
uses  the  equivalent  linear  algorithm  to  account  for  soil  nonlinearity.  The 
advantages  of  FLUSH  are  that  it  is  freely  downloadable  from  the  Internet  and  has 
considerably  faster  run  times  than  FLAG.  However,  the  major  disadvantage  of 
FLUSH  is  that  it  does  not  allow  for  permanent  displacement  of  the  wall.  FLUSH 
accounts  for  the  nonlinear  response  of  soils  during  earthquake  shaking  through 
adjustments  of  the  soil  (shear)  stiffhess  and  damping  parameters  (as  a  function  of 
shear  strain)  that  develop  in  each  element  of  the  finite  element  mesh.  The  FLAG 
and  FLUSH  results  will  be  compared. 
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2  Selection  of  Design  Ground 
Motion 


2.1  Selection  Criteria 

The  selection  of  an  earthquake  acceleration  time-histoiy  for  use  in  the 
numerical  analyses  was  guided  by  the  following  criteria: 

a.  A  real  earthquake  motion  was  desired,  not  a  synthetic  motion. 

b.  The  earthquake  magnitude  and  site-to-source  distance  corresponding  to 
the  motion  should  be  representative  of  design  ground  motions. 

c.  The  motion  should  have  been  recorded  on  rock  or  stiff  soil. 

These  criteria  were  used  to  assemble  a  list  of  candidate  acceleration  time- 
histories,  while  additional  criteria,  discussed  in  Section  2.3,  were  used  to  select 
one  time-history  from  the  candidate  list.  Because  the  response  of  a  soil-structure 
system  in  a  linear  dynamic  analysis  is  governed  primarily  by  the  spectral  content 
of  the  time-history  and  because  it  is  possible  to  obtain  a  very  close  fit  to  the 
design  spectrum  using  spectrum-matching  methods,  it  is  sufficient  to  have  a 
single  time-history  for  each  component  of  motion  for  each  design  earthquake. 
However,  because  the  nonlinear  response  of  a  soil-structure  system  may  be 
strongly  affected  by  the  time-domain  character  of  the  time-histories  even  if  the 
spectra  of  different  time-histories  are  nearly  identical,  at  least  five  time-histories 
(for  each  component  of  motion)  should  be  used  for  each  design  earthquake 
(Engineering  Circular  (EC)  1110-2-6051  (HQUSACE  2000)).  More  time- 
histories  are  required  for  nonlinear  dynamic  analyses  than  for  linear  analyses 
because  the  dynamic  response  of  a  nonlinear  structure  may  be  importantly  influ¬ 
enced  by  the  time  domain  character  of  the  time-history  (e.g.,  shape,  sequence, 
and  number  of  pulses),  in  addition  to  the  response  spectrum  characteristics. 
However,  for  the  first  phase  of  this  research  investigation,  only  one  time-history 
was  selected  for  use  in  the  dynamic  analyses. 


2.1.1  Real  versus  synthetic  earthquake  motion 

Because  the  numerical  analyses  performed  in  the  first  phase  of  this  research 
investigation  involve  permanent  displacement  of  the  wall  and  plastic  deforma¬ 
tions  in  the  soil  (i.e.,  nonlinearity),  it  was  decided  that  a  real  motion  should  be 
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used.  The  rationale  for  this  decision  was  to  avoid  potential  problems  of  develop¬ 
ing  a  synthetic  motion  that  appropriately  incorporates  all  the  factors  that  may 
influence  the  dynamic  response  of  a  nonlinear  system. 


2.1.2  Representative  magnitude  and  site-to-source  distance 

As  stated  in  Chapter  1,  the  objective  of  this  study  is  to  determine  the  seismic 
structural  design  loads  for  the  stem  portion  of  a  cantilever  retaining  wall. 
Accordingly,  the  magnitude  Mand  site-to-source  distance  R  of  the  ground 
motion  used  in  the  numerical  analyses  should  be  representative  of  an  actual 
design  earthquake,  which  will  depend  on  several  factors  including  geographic 
location  and  consequences  of  failure.  In  an  effort  to  select  a  "representative"  M 
and  R  for  a  design  event,  the  deaggregated  hazard  of  five  cities  located  in  the 
western  United  States  (WUS)  were  examined:  San  Francisco,  Oakland, 

Los  Angeles,  San  Diego,  and  Salt  Lake  City.  Deaggregation  of  the  seismic 
hazard  is  a  technique  used  in  conjunction  with  probabilistic  seismic  hazard 
analyses  (PSHA)  (EM  1 1 10-2-6050  (HQUSACE  1999))  to  express  the 
contribution  of  various  A/ and  R  combinations  to  the  overall  seismic  hazard  at  a 
site.  The  deaggregation  results  are  often  described  in  terms  of  the  mean 

magnitude  M  and  mean  distance  R  for  various  spectral  frequencies  (Frankel  et  al. 
1997).  It  is  not  uncommon  to  set  the  design  earthquake  magnitude  and  distance 

equal  to  the  values  of  M  and  R  corresponding  to  the  fundamental  frequency  of 
the  system  being  designed. 


Table  2-1  lists  the  M  and  R  for  the  peak  ground  acceleration  pga  and  1-hz 
spectral  acceleration  for  the  five  WUS  cities.  These  ground  motions  have  aver¬ 
age  return  periods  of  about  2500  years  (i.e.,  2  percent  probability  of  exceedance 
in  50  years).  From  the  deaggregated  hazards,  representative  Mand  R  for  the 
design  ground  motions  were  selected  as  7.0  and  25  km,  respectively. 


Table  2-1 

Mean  Magnitudes  a 
2500-year  Ground  1 

ind  Distances  for  Five  WUS  Cities  for  the 

I/lotion 

WUS  City 

pga 

Rpga .  km 

^ihz  -  km 

San  Francisco,  CA 

7.8 

25.0 

7.9 

25.0 

Oakland,  CA 

7.2 

25.0 

7.3 

25.4 

Los  Angeles,  CA 

6.8 

25.2 

7.0 

27.1 

San  Diego,  CA 

7.0 

25.0 

7.0 

25.1 

Salt  Lake  City.  UT 

7.1 

25.1 

7.3 

25.1 

2.1.3  Site  characteristics  of  motion 

The  amplitude  and  frequency  content,  as  well  as  the  phasing  of  the  frequen¬ 
cies,  of  recorded  earthquake  motions  are  influenced  by  the  source  mechanism 
(i.e.,  fault  type  and  rupture  process),  travel  path,  and  local  site  conditions,  among 
other  factors.  Because  the  selected  ground  motion  ultimately  is  to  be  specified  as 
a  base  rock  motion  in  the  numerical  analyses,  the  site  condition  for  the  selected 
ground  motions  is  desired  to  be  as  close  as  possible  to  the  base  rock  conditions 
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underlying  the  profile  on  which  the  cantilever  wall  is  located.  This  avoids  addi¬ 
tional  processing  of  the  recorded  motion  to  remove  the  site  effects  on  which  it 
was  recorded  (e.g.,  deconvolving  the  record  to  base  rock).  Accordingly,  motions 
recorded  on  rock  or  stiff  soil  profiles  were  desired  for  this  study. 


2.2  List  of  Candidate  Motions 


Based  on  the  selection  criteria,  the  motions  listed  in  Table  2-2  were 
considered  as  candidates  for  use  in  the  numerical  analyses. 


Table  2-2 
Candidate  A 

/lotions 

Earthquake 

Station 

Record 

pga.g 

Cape 

89530  Shelter  Cove  Airport 

SHL-UP 

0.054 

Mendocino 

Closest  to  fault  rupture:  33.8  km 

SHLOOO 

0.229 

M7.1,  Ms7.1 

Closest  to  surface  projection  of  rupture:  32.6  km 

SHL090 

0.189 

Duzce,  Turkey 

1058  Lamont 

1058-E 

0.111 

M7.1,  Ms7.3 

Closest  to  fault  rupture:  0.9  km 

1058-N 

0.073 

Closest  to  surface  projection  of  rupture:  0.9  km 

1058-V 

0.07 

Duzce,  Turkey 

1061  Lamont 

1061-E 

0.134 

M7.1,  Ms7.3 

Closest  to  fault  rupture:  15.6  km 

1061-N 

0.107 

Closest  to  surface  projection  of  rupture:  15. 6  km 

1061-V 

0.048 

Loma  Prieta 

57383  Gilroy  Array  #6 

G06-UP 

0.101 

M6.9,  Ms7.1 

Closest  to  fault  rupture:  19.9  km 

G06000 

0.126 

Closest  to  surface  projection  of  rupture:  19.9  km 

G06000 

0.1 

Loma  Prieta 

47189  SAGO  South-surface 

SG3-UP 

0.06 

M6.9,  Ms7.1 

Closest  to  fault  rupture:  34. 7  km 

SG3261 

0.073 

Closest  to  surface  projection  of  rupture:  34. 1  km 

SG3351 

0.067 

1  Note:  Ms  =  surface  wave  magnitude  of  earthquake;  M  =  moment  mai 

gnitude  of  earthquake.  I 

These  records  were  obtained  by  searching  the  Strong  Motion  Database 
maintained  by  the  Pacific  Earthquake  Engineering  Research  (PEER)  Center 
(http://peer.berkelev.edu/smcat/). 


2.3  Characteristics  of  Ground  Motion  Selected 

As  stated  previously,  at  least  five  time-histories  (for  each  component  of 
motion)  meeting  the  selection  criteria  should  be  used  in  nonlinear  dynamic  analy¬ 
ses  (EC  1 1 10-2-6051  (HQUSACE  2000)).  However,  for  the  first  phase  of  this 
study,  only  SG3351  was  used,  which  was  recorded  during  the  1989  Loma  Prieta 
earthquake  in  California.  The  basis  for  selecting  SG335 1  was  that  it  was  esti¬ 
mated,  using  CWROTATE  (Ebeling  and  White,  in  preparation),  to  induce  the 
greatest  permanent  relative  displacement  of  the  wall.  The  numerical  formulation 
in  CWROTATE  is  based  on  the  Newmark  sliding  block  procedure  outlined  in 
Ebeling  and  Morrison  (1992),  Section  6.3,  and  is  discussed  further  in 
Appendix  C. 

SG3351  is  plotted  in  Figure  2-1,  as  well  as  the  corresponding  5  percent 
damped,  pseudo-acceleration  response  spectrum,  scaled  to  1  g  pga.  Additionally, 
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Pseudo  Spectral  Acceleration  /  pga 


Figure  2-2.  Husid  plot  of  SG3351  used  for  determining  duration  of  strong  shaking, 
18.3  sec  (a(Q  is  the  acceleration  at  time  t  and  tf  is  the  total  duration  of 


the  acceleration  time-history) 


a  Husid  plot  of  the  motion  is  shown  in  Figure  2-2,  which  was  used  to  compute 
duration  of  strong  shaking  (EC  1 1 10-2-605 1  (HQUSACE  2000)),  1 8.3  sec. 


2.4  Processing  of  the  Selected  Ground  Motion 

Although  motion  SG335 1  met  the  selection  criteria,  several  stages  of 
processing  were  required  before  it  could  be  used  as  an  input  motion  in  the  FLAG 
analyses.  The  first  stage  was  simply  scaling  the  record.  As  a  general  rule, 
ground  motions  can  be  scaled  upward  by  a  factor  of  two  without  distorting  the 
realistic  characteristics  of  the  motion  (EC  1 1 10-2-6051  (HQUSACE  2000)).  The 
upward  scaling  was  desired  because  although  the  motion  induced  the  largest 
permanent  relative  displacement  dr  of  the  candidate  records,  the  induced 
displacement  was  still  too  small  to  ensure  active  earth  pressures  were  achieved. 
For  the  retaining  wall  system  being  modeled  in  this  first  phase  (i.e.,  wall  height: 
20  ft  (6  m);  backfill:  medium-dense,  compacted)  dr  >  0.5  in.  (12.7  mm)  is 
required  for  active  earth  pressures  (Ebeling  and  Morrison,  1992,  Table  1,  as 
adapted  from  values  presented  in  Clough  and  Duncan  1991). 

The  second  processing  stage  involved  filtering  high  frequencies  and 
computing  the  corresponding  interlayer  motion,  both  of  which  are  required  for 
either  finite  element  or  finite  difference  analyses.  As  discussed  subsequently  in 
detail  in  Chapter  3,  in  the  finite  element  and  finite  difference  formulations,  the 
element  dimension  in  the  direction  of  wave  propagation,  as  well  as  the 
propagation  velocity  of  the  material,  limits  the  maximum  frequency  which  the 
element  can  accurately  transfer.  For  most  soil  systems  and  most  earthquake 
motions,  the  removal  of  frequencies  above  15  hz  (i.e.,  low-pass  cutoff  frequency) 
will  not  influence  the  dynamic  response  of  the  system.  However,  caution  should 
be  used  in  selecting  the  low-pass  cutoff  frequency,  especially  when  the  motions 
are  being  used  in  dynamic  soil-structure-interaction  analyses  where  the  building 
structure  may  have  a  high  natural  frequency,  such  jis  nuclear  power  plants.  Next, 
SG3351  was  recorded  on  the  ground  surface,  and  the  corresponding  interlayer 
motion  needed  to  be  computed  for  input  into  the  base  of  the  FLAG  model.  A 
modified  version  of  the  computer  program  SHAKE91  (Idriss  and  Sun  1992)  was 
used  both  to  remove  the  high  frequencies  and  compute  the  interlayer  motion. 
Figure  2-3  shows  the  recorded  SG335 1  and  the  processed  record  used  as  input  at 
the  base  of  the  FLAG  model. 
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Figure  2-3.  Selected  ground  motion  (a)  recorded  motion  SG3351and  (b)  the 
processed  motion  used  as  input  into  the  base  of  the  FLAG  model 
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3  Numerical  Analysis  of 
Cantilever  Retaining  Wall 


3.1  Overview  of  FLAG 

As  stated  in  Chapter  1,  the  detailed  numerical  analyses  of  the  cantilever 
retaining  walls  were  performed  using  FLAG,  a  commercially  available,  two- 
dimensional,  explicit  finite  difference  program,  which  was  written  primarily  for 
geotechnical  engineering  applications.  The  basic  formulation  of  FLAG  is  plane- 
strain,  which  is  the  condition  associated  with  long  structures  perpendicular  to  the 
analysis  plane  (e.g.,  retaining  wall  systems).  The  following  is  a  brief  overview  of 
FLAG  and  is  largely  based  on  information  provided  in  the  FLAG  manuals  (Itasca 
Consulting  Group,  Inc.,  2000).  The  reader  is  referred  to  these  manuals  for 
additional  details. 

Because  it  is  likely  that  most  readers  are  more  familiar  with  the  finite  element 
method  (FEM)  than  with  the  finite  difference  method  (FDM),  analogous  terms  of 
the  two  methods  are  compared  as  shown: 

Finite  Element  Finite  Difference 


element 

zone 

node 

O 

grid  point 

mesh 

<=> 

grid 

In  places  of  convenience,  these  terms  are  used  interchangeably  in  this  report. 
For  example,  the  terms  structural  elements  and  interface  elements  are  used  in 
this  report,  as  opposed  to  structural  zones  and  interface  zones.  Both  FEM  and 
FDM  translate  a  set  of  differential  equations  into  matrix  equations  for  each 
element,  relating  forces  at  nodes  to  displacements  at  nodes.  The  primary 
difference  between  FLAG  and  most  finite  element  programs  relates  to  the 
explicit,  Lagrangian  calculation  scheme  used  in  FLAG,  rather  than  the 
differences  between  the  FEM  and  FDM.  However,  neither  the  Lagrangian 
formulation  nor  the  explicit  solution  scheme  is  inherently  unique  to  the  FDM  and 
may  be  used  in  the  FEM. 
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Dynamic  analyses  can  be  performed  with  FLAG  using  the  optional  dynamic 
calculation  module,  wherein  user-specified  acceleration,  velocity,  or  stress  time- 
histories  can  be  input  as  an  exterior  boundary  condition  or  as  an  interior 
excitation.  FLAG  allows  energy-absorbing  boundary  conditions  to  be  specified, 
which  limits  the  numerical  refiection  of  seismic  waves  at  the  model  perimeter. 

FLAG  has  ten  built-in  constitutive  models,  including  a  null  model,  and 
allows  user-defined  models  to  be  incorporated.  The  null  model  is  commonly 
used  in  simulating  excavations  or  construction,  where  the  finite  difference  zones 
are  assigned  no  mechanical  properties  for  a  portion  of  the  analysis.  The  explicit 
solution  scheme  can  follow  arbitrary  nonlinear  stress-strain  laws  with  little 
additional  computational  effort  over  linear  stress-strain  laws.  FLAG  solves  the 
full  dynamic  equations  of  motion,  even  for  essentially  static  systems,  which 
enables  accurate  modeling  of  unstable  processes  (e.g.,  retaining  wall  failures). 
The  explicit  calculation  cycle  used  in  FLAG  is  illustrated  in  Figure  3-1. 


Figure  3-1 .  Basic  explicit  calculation  cycle  used  in  FLAG  (adapted  from  Itasca 
Consulting  Group,  Inc.,  2000,  Theory  and  Background  Manual) 


Referring  to  Figure  3-1,  beginning  with  a  known  stress  state,  the  equation  of 
motion  is  solved  for  the  velocities  and  displacements  for  each  element,  while  it  is 
assumed  that  the  stresses  are  frozen.  Next,  using  the  newly  computed  velocities 
and  displacements,  in  conjunction  with  the  specified  stress-strain  law,  the  stresses 
are  computed  for  each  element,  while  it  is  assumed  that  the  velocities  and 
displacements  are  frozen.  The  assumption  of  frozen  velocities  and  displacements 
while  stresses  are  computed  and  vice-versa  can  produce  accurate  results  only  if 
the  computational  cycle  is  performed  for  a  very  small  increment  in  time  (i.e.,  the 
"calculation  wave  speed"  must  always  be  faster  than  the  physical  wave  speed). 
This  leads  to  the  greatest  disadvantage  of  FLAG,  long  computational  times, 
particularly  when  modeling  stiff  materials,  which  have  large  physical  wave 
speeds.  The  size  of  the  time-step  depends  on  the  dimension  of  the  elements,  the 
wave  speed  of  the  material,  and  the  type  of  damping  specified  (i.e.,  mass 
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proportional  or  stiffness  proportional),  where  stiffness  proportional,  to  include 
Rayleigh  damping,  requires  a  much  smaller  time-step.  The  critical  time-step  for 
stability  and  accuracy  considerations  is  automatically  computed  by  FLAG,  based 
on  these  factors  listed.  For  those  readers  unfamiliar  with  the  concept  of  critical 
time-step  for  stability  and  accuracy  considerations  in  a  seismic  time-history 
engineering  analysis  procedure,  please  refer  to  Ebeling  (1992),  Part  V,  or  to 
Ebeling,  Green,  and  French  (1997). 

The  Lagrangian  formulation  in  FLAG  updates  the  grid  coordinates  each 
time-step,  thus  allowing  large  cumulative  deformations  to  be  modeled.  This  is  in 
contrast  to  the  Eulerian  formulation  in  which  the  material  moves  and  deforms 
relative  to  a  fixed  grid,  and  is  therefore  limited  to  small  deformation  analyses. 


3.2  Retaining  Waii  Modei 

The  retaining  wall-soil  system  analyzed  in  the  first  phase  of  this  investigation 
is  depicted  in  Figure  3-2.  As  shown  in  this  figure,  the  FLAG  model  is  only  the 
top  30  ft  (9  m)  of  a  225-ft  (69-m)  profile.  Although  the  entire  profile,  to  include 
the  retaining  wall,  can  be  modeled  in  FLAG,  the  required  computational  time 
would  be  exorbitant,  with  little  to  no  benefit  added.  To  account  for  the  influence 
of  the  soil  profile  below  30  ft  (9  m),  the  entire  profile  without  the  retaining  wall 
was  modeled  using  a  modified  version  of  SHAKE91  (Idriss  and  Sun  1992),  and 
the  interlayer  motion  at  the  depth  corresponding  to  the  base  of  the  FLAG  model 
was  computed.  The  input  ground  motion  used  in  the  SHAKE  analysis  was 
SG3351,  the  selection  of  which  was  discussed  in  Ghapter  2.  SG3351  was 
specified  as  a  rock  outcrop  motion  for  the  soil  column.  In  this  type  of  analysis 
the  base  of  the  soil  column  is  modeled  as  a  halfspace  in  the  SHAKE  model.  In 
order  to  account  for  the  site-specific  pga  value  anticipated  at  this  site  for  the 
specified  design  earthquake  magnitude  and  specified  site-to-source  distance 
(discussed  in  Ghapter  2),  a  scale  factor  of  two  was  applied  to  SG3351 
acceleration  time-history.  Based  on  the  guidelines  in  EG  1 1 10-2-6051 
(HQUSAGE  2000)  allowing  motions  to  be  scaled  upward  by  a  factor  less  than  or 
equal  to  two,  this  action  was  judged  to  be  reasonable  by  this  Gorps  criterion.  The 
variation  of  the  shear  wave  velocity  as  a  function  of  depth  in  the  profile  is 
consistent  with  dense  natural  deposits  beneath  the  base  of  the  retaining  wall  and 
medium-dense  compacted  fill  for  the  backfill. 

The  small  strain  natural  frequency  of  the  retaining  wall-soil  system  in  the 
FLAG  model  is  estimated  to  be  approximately  6.2  hz,  as  determined  by  the  peak 
of  the  transfer  function  from  the  base  of  the  model  to  the  top  of  the  backfill.  At 
higher  strains,  it  is  expected  that  the  natural  period  of  the  system  will  be  less  than 

6.2  hz.  The  cutoff  frequency  in  the  SHAKE  analysis  was  set  at  15  hz.  This 
value  was  selected  based  on  both  the  natural  frequency  of  the  wall-soil  system 
and  the  energies  associated  with  the  various  frequencies  in  SG3351,  and  ensures 
proper  excitation  of  the  wall.  Dimensioning  of  the  finite  difference  zones  to 
ensure  proper  transfer  of  frequencies  up  to  15  hz  is  discussed  in  Section  3.3.4. 
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Figure  3-2.  Numerical  models  used  in  the  dynamic  analysis  of  the  cantilever  retaining  wall 
(To  convert  feet  to  meters,  multiply  by  0.3048) 
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The  interlayer  motion  computed  using  SHAKE  was  specified  as  an  accelera¬ 
tion  time-history  along  the  base  of  the  FLAG  model.  Based  on  the  results  of  a 
parametric  study,  outlined  in  Appendix  D,  specification  of  the  ground  motion  in 
FLAG  in  terms  of  acceleration,  as  opposed  to  stress  or  velocity,  gives  the  most 
accurate  results  for  the  profiles  analyzed. 

Figure  3-3  shows  an  enlargement  of  the  retaining  wall-soil  system  modeled 
in  FLAG,  as  well  as  the  finite  difference  grid.  The  FLAG  model  consists  of  four 
subgrids,  labeled  1  through  4.  Subgrids  are  used  in  FLAG  to  create  regions  of 
different  shapes;  there  is  no  restriction  on  the  variation  of  the  material  properties 
of  the  zones  within  a  subgrid.  The  separation  of  the  foundation  soil  and  backfill 
into  Subgrids  1  and  2  was  required  because  a  portion  of  the  base  of  the  retaining 
wall  is  inserted  into  the  soil.  Subgrid  4  was  required  because  the  free-field 
boundary  conditions,  an  energy-absorbing  boundaiy  option  in  FLAG,  cannot  be 
specified  across  the  interface  of  two  subgrids.  Subgrid  3  was  included  for 
symmetiy.  The  subgrids  were  “attached”  at  the  soil-to-soil  interfaces,  as  depicted 
by  the  dashed  red  line  in  Figure  3-3a,  and  the  yellow  +'s  in  Figure  3-3b.  The 
attach  command  welds  the  corresponding  grid  points  of  two  subgrids  together. 
Interface  elements  were  used  at  the  soil-structure  interfaces,  as  depicted  by  green 
lines  in  Figure  3-3a,  and  discussed  in  more  detail  in  Section  3.3.3. 


25  ft 


-*+*- 


50  ft 


b) 


Figure  3-3.  Retaining  wall-soil  system  modeled  in  FLAG:  (a)  conceptual  drawing 
showing  dimensions  and  soil  layering  and  (b)  finite  difference  grid 
(To  convert  feet  to  meters,  multiply  by  0.348) 
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The  retaining  wall  model  was  "numerically  constructed"  in  FLAG  similar  to 
the  way  an  actual  wall  would  be  constructed.  The  backfill  was  placed  in  2-ft 
(0.6-m)  lifts,  for  a  total  of  ten  lifts,  with  the  model  being  brought  to  static 
equilibrium  after  the  placement  of  each  lift.  This  allowed  realistic  earth  pressures 
to  develop  as  the  wall  deformed  and  moved  due  to  the  placement  of  each  lift. 
Figure  3-4  shows  the  deformed  grid,  magnified  75  times,  after  the  construction  of 
the  wall  and  backfill  placement  to  a  height  of  20  ft  (6  m). 

In  the  next  section,  an  overview  is  given  on  how  the  numerical  model 
parameters  were  determined. 


Figure  3-4.  Deformed  finite  difference  grid,  magnified  75  times.  The  backfill  was 
placed  in  ten  2-ft  (0.6-m)  lifts,  with  the  model  being  brought  to  static 
equilibrium  after  the  placement  of  each  lift 


3.3  Numerical  Model  Parameters 

The  previous  section  gave  an  overview  of  the  physical  system  being  analyzed 
and  its  numerical  model  counterpart.  This  section  focuses  on  the  specific 
constitutive  models  used  for  the  soil,  retaining  wall,  and  their  interface,  with 
particular  attention  given  to  how  to  determine  the  veirious  model  parameters.  An 
elastoplastic  constitutive  model,  in  conjunction  with  Mohr-Coulomb  failure 
criteria,  was  used  to  model  the  soil.  Elastic  beam  elements  were  used  to  model 
the  concrete  retaining  wall,  and  interface  elements  were  used  to  model  the 
interaction  between  the  soil  and  the  structure.  The  following  sections  outline  the 
procedures  used  to  determine  the  various  model  parameters. 


3.3.1  Mohr-Coulomb  model 

Four  parameters  are  required  for  the  Mohr-Coulomb  model:  internal  friction 
angle  <|);  mass  density  p;  shear  modulus  G;  and  bulk  modulus  K'.  The  first  two 
parameters,  (|)  and  p,  are  familiar  to  geotechnical  engineers,  where  mass  density  is 
the  total  unit  weight  of  the  soil  y,  divided  by  the  acceleration  due  to  gravity  g, 
i.e.,  p  =  y,/g.  (|)  for  the  foundation  soil  was  set  at  40  deg  and  to  35  deg  for  the 
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backfill.  These  values  are  consistent  with  dense  natural  deposits  and  medium- 
dense  compacted  fill.  G  and  K'  may  be  less  familiar  to  geotechnical  engineers, 
and  therefore,  their  determination  is  outlined  as  follows. 

Shear  modulus  G.  Several  correlations  exist  that  relate  G  to  other  soil 
parameters.  However,  the  most  direct  relation  is  between  G  and  shear  wave 
velocity  v^: 

G  =  pv]  (3-1) 

Vj  may  be  determined  by  various  types  of  site  characterization  techniques,  such  as 
cross  hole  or  spectral  analysis  of  surface  waves  (SASW)  studies. 

Bulk  modulus  K'.  Values  for  K'  are  typically  computed  from  G  and 
Poisson's  ratio  v  using  the  following  relation: 


2-G(l  +  v)  (3-2) 


For  natural  deposits,  v  may  be  estimated  from  the  following  expression: 


^  l-sin(<^)  (3-3) 

2-sin(^) 

This  expression  can  be  derived  from  the  theory  of  eleisticity  (Terzaghi  1943) 
and  the  correlation  relating  at-rest  earth  pressure  conditions  Ko  and  (|)  proposed  by 
Jaky  (1944).  However,  for  surficial  compacted  soil  against  a  nondeflecting  soil 
structure  interface,  Duncan  and  Seed  (1986)  proposed  the  following  ’’empirically 
derived”  expression  for  v,  which  results  in  considerably  higher  values  than  that 
for  natural  soil  deposits: 


4“3sin  (^) 
8-4*sin  (^) 


(3-4) 


For  the  numerical  analysis  of  a  retaining  wall  with  a  compacted  backfill,  for 
which  laboratory  tests  are  not  performed  to  determine  v,  judgment  should  be 
used  in  selecting  a  value  for  v,  with  a  reasonable  value  being  between  those 
given  by  Equations  3-3  and  3-4: 
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2-sin(^)  8-4-sin(^) 


(3-5) 


3.3.2  Structural  elements 

The  concrete  retaining  wall  was  modeled  using  elastic  beam  elements 
approximately  1  ft  (0.3  m)  long.  In  FLAG,  four  parameters  are  required  to  define 
the  mechanical  properties  of  the  beam  elements:  cross-sectional  area  Agi  mass 
density  p;  elastic  modulus  Ec',  and  second  moment  of  area  /,  commonly  referred 
to  as  moment  of  inertia.  The  wall  was  divided  into  five  segments  having  con¬ 
stant  parameters,  as  illustrated  in  Figure  3-5,  with  each  segment  consisting  of 
several  1-ft  (0.3-m)  beam  elements.  The  number  of  segments  used  is  a  function 
of  the  variation  of  the  mechanical  properties  in  the  wall.  A  wall  having  a  greater 
taper  or  largely  varying  reinforcing  steel  along  the  length  of  the  stem  or  base 
would  likely  require  more  segments. 


Figure  3-5.  Subdivision  of  the  cantilever  wall  into  five  segments,  each  having 
constant  material  properties  (To  convert  feet  to  meters,  multiply  by 
0.3048) 


For  each  of  the  segments,  Ag  and  p  were  readily  determined  from  the  wall 
geometry  and  the  unit  weight  of  the  concrete  (i.e.,  150  Ib/cu  ft  (2,403  kg/cu  m)). 
Ag  and  p  are  used  in  FLAG  to  compute  gravity  and  inertial  forces. 

Ec  was  computed  using  the  following  expression  (MacGregor  1992): 


Chapter  3  Numerical  Analysis  of  Cantilever  Retaining  Wall 


21 


(3-6) 


£,=57,000-^ 

In  this  expression,/,;  is  the  compressive  strength  of  the  concrete  (e.g., 

4,000  psi  (28  MPa)  for  the  wall  being  modeled),  and  both  Ec  and/c  are  in  psi. 
Because  the  structure  is  continuous  in  the  direction  perpendicular  to  the  analysis 
plane,  Ec  needs  to  be  modified  using  the  following  expression  to  account  for 
plane-strain  conditions,  where  0.2  was  assumed  for  Poisson's  ratio  for  concrete 
(Itasca  Consulting  Group,  Inc.,  2000,  FLAG  Structural  Elements  Manual): 


E. 


'c  plane  Strain 


Ec 


(3-7) 


I  is  dependant  on  the  geometry  of  the  segments,  the  amount  of  reinforcing 
steel,  and  the  amount  of  cracking  in  the  concrete,  where  the  latter  is  in  turn  a 
function  of  the  static  and  dynamic  loading  imposed  on  the  member.  Table  3-1 
presents  values  of  /for  the  five  wall  segments  assuming  uncracked  and  fully 
cracked  sectional  properties.  In  dynamic  analyses,  it  is  difficult  to  state  a  priori 
whether  use  of  sectional  properties  corresponding  to  uncracked,  fully  cracked,  or 
some  intermediate  level  of  cracking  will  result  in  the  largest  demand  on  the 
structure.  In  this  first  phase  of  this  research  investigation,  two  FLAG  analyses 
were  performed  assuming  the  extreme  values  for  /(i.e.,  hncracked  and  hracked)- 
However,  using  /=  (^A-Iuncracked  is  a  reasonable  estimate  for  the  sectional 
properties  for  most  cases  (Paulay  and  Priestley  1992). 


Table  3-1 

Second  Moment  of  Area  for  Wall  Segments 

Section 

funcracAtd  (ft  ) 

1 

0.339 

0.094 

2 

0.439 

0.121 

3 

0.557 

0.152 

4 

0.694 

0.186 

5 

0.701 

0.174 

Note:  To  convert  ft^  to  m\  multi( 

3ly  by  0.008631.  II 

3.3.3  Interface  elements 

Interface  elements  were  used  to  model  the  interaction  between  the  concrete 
retaining  wall  and  the  soil.  However,  FLAG  does  not  allow  interface  elements  to 
be  used  at  the  intersection  of  branching  structures  (e.g.,  the  intersection  of  the 
stem  and  base  of  the  cantilever  wall).  Of  the  several  attempts  by  the  authors  to 
circumvent  this  limitation  in  FLAG,  the  simplest  and  best  approach  found  is 
illustrated  in  Figure  3-6.  As  shown  in  this  figure,  three  very  short  beam  ele¬ 
ments,  oriented  in  the  direction  of  the  stem,  toe  side  of  the  base,  and  heel  side  of 
the  base,  were  used  to  model  the  base-stem  intersection.  No  interface  elements 
were  used  on  these  three  beam  elements.  However,  interface  elements  were  used 
along  the  other  contact  surfaces  between  the  soil  and  wall,  as  depicted  by  the 
green  lines  in  Figure  3-6. 
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h-  3ft  8ft  -►! 


Figure  3-6.  Approach  to  circumventing  the  limitation  in  FLAG  of  not  allowing 

interface  elements  to  be  used  at  branching  intersections  of  structural 
elements  (To  convert  feet  to  meters,  multiply  by  0.3048) 


A  schematic  of  the  FLAG  interface  element  and  the  inclusive  parameters  is 
presented  in  Figure  3-7.  The  element  allows  permanent  separation  and  slip  of  the 
soil  and  the  structure,  as  controlled  by  the  parameters  tensile  strength  T  and 
slider  S,  respectively.  For  the  analyses  performed  as  part  of  this  research 
investigation,  T  =  0,  thus  modeling  a  cohesionless  soil.  S  was  specified  as  a 
function  of  the  interface  friction  angle  6.  Based  on  the  values  of  8  for  medium- 
dense  sand  against  concrete  given  in  Tables  3-7  and  3-8  of  Gomez,  Filz,  and 
Ebeling  (2000b),  5  =  31  deg  was  selected  as  a  representative  value. 

Normal  stiffness  k„.  The  FLAG  manual  (Itasca  Gonsulting  Group,  Inc., 
2000,  Theory  and  Background  Manual)  recommends  as  a  rule  of  thumb  that  be 
set  to  ten  times  the  equivalent  stiffness  of  the  stiffest  neighboring  zone,  i.e.. 


k 

n 


10 -max 


4 

K'  +  —  G 
3 


Az  . 
min 


(3-8) 


In  this  relation,  Az„,„  is  the  smallest  width  of  a  zone  in  the  normal  direction 
of  the  interfacing  surface.  The  max[  ]  notation  indicates  that  the  maximum  value 
over  all  zones  adjacent  to  the  interface  is  used.  The  FLAG  manual  warns  against 
using  arbitrarily  large  values  for  k„,  as  is  commonly  done  in  finite  element 
analyses,  as  this  results  in  an  unnecessarily  small  time-step  and  therefore 
unnecessarily  long  computational  times. 
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Figure  3-7.  Schematic  of  the  FLAG  interface  element  (adapted  from  Itasca 
Consulting  Group,  Inc.,  2000,  Theory  and  Background  Manual) 


Shear  stiffness  k,.  The  determination  of  ks  required  considerably  more  effort 
than  the  determination  of  the  other  interface  element  parameters.  In  shear,  the 
interface  element  in  FLAG  essentially  is  an  elastoplastic  model  with  an  elastic 
stiffness  of  ks  and  yield  strength  S.  ks  values  were  selected  such  that  the  resulting 
elasto-plastic  model  gave  an  approximate  fit  of  the  hyperbolic-type  interface 
model  proposed  by  Gomez,  Filz,  and  Ebeling  (2000a,b).  A  comparison  of  the 
two  models  is  shown  in  Figure  3-8  for  initial  loading  (i.e.,  construction  of  the 
wall). 

The  procedure  used  to  determine  ks  values  for  initial  loading  is  outlined  in 
the  following  steps.  See  Gomez,  Filz,  and  Ebeling  (2000a,b)  for  more  details 
concerning  their  proposed  hyperbolic-type  model. 

a.  Compute  the  reference  displacement  along  the  interface  Ar  using  the 
following  expression.  Representative  values  for  Rg,  Ki,  rij,  and  5  were 
obtained  from  Gomez,  Filz,  and  Ebeling  (2000a). 


Ar  = 


(3-9a) 


where 


=cr„-tan(J) 


(3-9b) 


(3-9c) 
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Figure  3-8.  Comparison  of  the  Gomez,  Filz,  and  Ebeling  (2000a, b)  hyperbolic-type  interface 

element  model  and  the  approximate-fit  elastoplastic  model  (To  convert  pounds  (force) 
per  square  foot  to  pascals,  multiply  by  47.88;  to  convert  feet  to  meters,  multiply  by 
0.3048) 

Rg  =  failure  ratio  =  0.84 

Ksi  =  initial  shear  stiffness  of  the  interface 

an  =  normal  stress  acting  on  the  interface,  and  determined  iteratively  in 
FLAG  by  first  assuming  a  small  value  for  kg  and  then  constructing  the 
wall 

6  =  interface  friction  angle  =  3 1  deg 

K[  =  dimensionless  interface  stiffness  number  for  initial  loading  =  21000 

y„,  =  unit  weight  of  water  in  consistent  units  as  Ar  (i.e.,  =  62.4  Ib/cu  ft 
(1,000  kg/cu  m)) 

Pa  =  atmospheric  pressure  in  the  same  units  as  a„ 
rij  =  dimensionless  stiffness  exponent  =  0.8 
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b.  ks  is  computed  using  the  following  expression: 


k,  =■ 


'a  J 


(3-10) 


These  computed  k^  values  were  used  only  for  the  initial  construction  of  the 
wall.  After  the  construction  of  the  wall  and  prior  to  the  application  of  the 
earthquake  loading  the  ks  values  were  changed  to  values  consistent  with  the 
Gomez-Filz-Ebeling  Version  I  load/unload/reload  extended  hyperbolic  interface 
model  (Gomez,  Filz,  and  Ebeling  2000b).  The  procedure  used  to  compute  ks  for 
the  cyclic  loading  is  outlined  in  the  following  equations.  Again,  refer  to  Gomez, 
Filz,  and  Ebeling  for  more  details  concerning  this  model. 


^s^^un-Yy. 


\^a  J 


(3-1  la) 


where 


K„=C,K,  <3-1  lb) 

C,  =0.5-(l  +  /!^)’  <3-1  li) 


and 


Kurj  =  unload-reload  stiffness  number  for  interfaces 


Ck  =  interface  stiffness  ratio 

The  interface  stiffnesses  were  computed  using  Equations  3-8,  3-10,  and  3-1 1 
for  the  interface  elements  identified  in  Figure  3-9.  The  computed  values  are 
listed  in  Table  3-2.  The  normal  stiffnesses  were  the  same  for  both  the  initial 
loading  and  the  unload-reload.  However,  the  shear  stiffnesses  increased  from  the 
initial  loading  to  the  unload-reload. 


3.3.4  Dimensions  of  finite  difference  zones 

As  mentioned  previously,  proper  dimensioning  of  the  finite  difference  zones 
is  required  to  avoid  numerical  distortion  of  propagating  ground  motions,  in  addi¬ 
tion  to  accurate  computation  of  model  response.  The  FLAG  manual  (Itasca 
Consulting  Group,  Inc.,  2000,  Optional  Features  Manual)  recommends  that  the 
length  of  the  element  A/  be  smaller  than  one-tenth  to  one-eighth  of  the  wave¬ 
length  \  associated  with  the  highest  frequency  fmax  component  of  the  input 
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Figure  3-9.  Interface  element  numbering 


Table  3-2 

k„  and  kg  Initial  and  Unload-reload  Values  for  Interface  Elements 

Interface  Element 

kn 

ks  initiaUoading 

ks  unload-reload 

Number 

Ibf/fl^/ft 

Ibf/ft^/ft 

ibf/ft'm 

1 

1.203x10“ 

7.875x10“ 

2.666x10“ 

2 

1.203x10“ 

6.709x10“ 

2.379x10“ 

3 

1.073x10“ 

9.967x10“ 

3.374x10“ 

4 

1.073x10“ 

9.012x10“ 

3.051x10“ 

5 

1.073x10“ 

3.056x10^ 

1.035x10“ 

6 

1.073x10“ 

2.647x10“ 

8.963x10“ 

7 

9.657x10' 

2.255x10“ 

7.634x10“ 

8 

9.657x10' 

1.827x10“ 

6.185x10“ 

9 

8.356x10' 

1.491x10“ 

5.050x10“ 

10 

8.356x10' 

1.214x10“ 

4.110x10“ 

11 

8.356x10' 

7.971x10" 

2.699x10“ 

12 

6.398x10' 

4.641x10* 

1.571x10“ 

13 

6.398x10' 

2.465x10* 

8.345x10* 

1  Note:  To  convert  Ibf/ft^/ft  to  pascals  per  meter,  multiply  by  157.09  || 

motion.  The  basis  for  this  recommendation  is  a  study  by  Kuhlemeyer  and 
Lysmer  (1973).  Interestingly,  the  FLUSH  manual  (Lysmer  et  al.  1975) 
recommends  A/  be  smaller  than  one-fifth  the  A,  associated  with  f„uix,  also 
referencing  Kuhlemeyer  and  Lysmer  (1973)  as  the  basis  for  the  recommendation, 
i.e.,: 

FMC  FLUSH 

A/<:^  A/<|  (3-12) 
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X  is  related  to  the  shear  wave  velocity  of  the  soil  and  the  frequency  /of  the 
propagating  wave  by  the  following  relation: 


A 


(3-13) 


In  a  FLUSH  analysis  it  is  important  to  note  that  the  used  in  this  compu¬ 
tation  is  not  that  for  small  (shear)  strains,  such  as  measured  in  the  field  using 
cross-hole  shear  wave  tests.  In  FLUSH,  the  used  to  dimension  the  elements 
should  be  consistent  with  the  earthquake-induced  shear  strains,  frequently 
referred  to  as  the  "reduced"  by  FLUSH  users.  Assuming  that  the  response  of 
the  retaining  wall  will  be  dominated  by  shear  waves,  substitution  of  Equa¬ 
tion  3-13  into  the  FLAG  expression  for  A/  in  Equation  3-12  gives: 


A/<— ^ 

io/„ 


(3-14a) 


f  ^ 

J  max  — 


10  Al 


(3-14b) 


As  may  be  observed  from  these  expressions,  the  finite  difference  zone  with 
the  lowest  and  a  given  Al  will  limit  the  highest  frequency  that  can  pass  through 
the  zone  without  numerical  distortion.  For  the  FLAG  analyses  performed  in 
Phase  1  of  this  investigation,  1-ft  by  1-ft  (0.3-m  by  0.3-m)  zones  were  used  in 
subgrids  1  and  2  (Figure  3-3).  The  top  layer  of  the  backfill  has  the  lowest  (i.e., 
525  fps  (160  m/sec)).  Using  these  expressions  and  A/  =  1  ft  (0.3  m),  the  finite 
difference  grid  used  in  the  FLAG  analyses  should  adequately  propagate  shear 
waves  having  frequencies  up  to  52.5  hz.  This  value  is  well  above  the  15-hz 
cutoff  frequency  used  in  the  SHAKE  analysis  to  compute  the  input  motion  for 
the  FLAG  analysis  and  well  above  the  estimated  fundamental  frequency  of  the 
retaining  wall-soil  system  being  modeled. 


3.3.5  Damping 

An  elastoplastic  constitutive  model  in  conjunction  with  the  Mohr-Goulomb 
failure  criterion  was  used  to  model  the  soil  in  FLAG.  Inherent  in  this  model  is 
the  characteristic  that  once  the  induced  dynamic  shear  stresses  exceed  the  shear 
strength  of  the  soil,  the  plastic  deformation  of  the  soil  introduces  considerable 
hysteretic  damping.  However,  for  dynamic  shear  stresses  less  than  the  shear 
strength  of  the  soil,  the  soil  behaves  elastically  (i.e.,  no  damping),  unless 
additional  mechanical  damping  is  specified.  FLAG  allows  mass  proportional, 
stiffness  proportional,  and  Rayleigh  damping  to  be  specified,  where  the  latter 
provides  relatively  constant  level  of  damping  over  a  restricted  range  of 
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frequencies.  Use  of  either  stiffness  proportional  or  Rayleigh  damping  results  in 
considerably  longer  run  times  than  when  either  no  damping  or  mass  proportional 
damping  is  specified. 

For  the  FLAG  analyses  performed,  Rayleigh  damping  was  specified,  for 
which  the  critical  damping  ratio  ^  may  be  determined  by  the  following  relation: 


a 


\G> 


+  P-O) 


(3-15) 


where 

a  =  the  mass-proportional  deimping  constant 
P  =  the  stiffness-proportional  damping  constant 
(0  =  singular  frequency  associated  with  ^ 

For  Rayleigh  damping,  the  damping  ratio  and  the  corresponding  central 
frequency  need  to  be  specified.  Judgment  is  required  in  selecting  values  for  both 
parameters.  A  lower  bound  for  the  damping  ratio  is  1  to  2  percent.  This  level 
helps  reduce  high-frequency  spurious  noise.  However,  considerable  high- 
frequency  noise  will  still  exist  even  when  1  to  2  percent  Rayleigh  damping  is 
specified;  this  is  an  inherent  shortcoming  of  an  explicit  solution  algorithm.  The 
damping  levels  in  the  last  iteration  of  SHAKE  analysis  used  to  compute  the 
FLAG  input  motion  may  be  used  as  an  upper  bound  of  the  values  for  Rayleigh 
damping.  The  central  frequency  corresponding  to  the  specified  damping  ratio  is 
typically  set  to  either  the  funcUimental  period  of  the  system  being  modeled  (an 
inherent  property  of  the  soil-wall  system)  or  predominant  period  of  the  system 
response  (an  inherent  property  of  the  soil-wall  system  and  the  ground  motion). 
For  the  FLAG  analyses  performed  as  part  of  the  first  phase  of  this  investigation, 
the  SHAKE  computed  damping  ratios  were  used  and  the  central  frequency  was 
set  equal  to  the  fundamental  frequency  of  the  retaining  wall-soil  system.  In 
future  FLAG  analyses,  the  damping  ratios  will  be  set  to  half  the  values  from  the 
SHAKE  analysis. 


3.4  Summary 

When  a  physical  system  is  modeled  numerically,  considerable  judgment  is 
required  in  selecting  appropriate  values  for  the  model  parameters.  This  chapter 
provided  an  overview  of  the  models  used  in  the  FLAG  analyses  of  the  cantilever 
retaining  wall  and  discussed  approaches  for  selecting  values  for  the  various 
model  parameters.  In  the  next  chapter,  the  results  from  the  FLAG  analyses  are 
presented  and  put  into  perspective  of  the  current  Gorps  design  procedure,  as 
presented  in  Ebeling  and  Morrison  (1992). 
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4  FLAG  Data  Reduction 
Discussion  of  Results 


In  the  previous  chapter,  an  overview  was  given  of  the  numerical  model  used 
to  analyze  the  cantilever  retaining  wall.  In  this  chapter  an  overview  is  given  on 
how  the  FLAG  data  were  reduced,  followed  by  a  presentation  and  discussion  of 
the  reduced  data.  Two  FLAG  analyses  were  performed  as  part  of  the  first  phase 
of  this  research  effort,  one  using  the  uncracked  properties  of  the  concrete  wall, 
and  the  other  using  the  fully  cracked  properties  (refer  to  Table  3-1  for  the  listing 
of  the  respective  properties).  The  results  from  the  two  analyses  were  similar. 
However,  more  information  was  computed  in  the  FLAG  analysis  of  the  fully 
cracked  wall  (i.e.,  additional  acceleration  time-histories  in  the  foundation  and 
backfill  soil  and  displacement  time-histories  along  the  wall  were  requested  in  the 
FLAG  input  file  for  the  cracked  wall).  Accordingly,  the  results  of  the  FLAG 
analysis  of  the  cracked  wall  are  primarily  discussed  in  this  chapter. 


4.1  Data  Reduction 

Time-histories  for  the  lateral  stresses  acting  on  the  elements  composing  the 
stem  and  heel  section  were  computed  by  FLAG,  as  well  as  acceleration  time- 
histories  at  various  places  in  the  foundation  soil  and  backfill  and  displacement 
time-histories  at  various  places  on  the  wall.  The  stresses  computed  by  FLAG  are 
averages  of  the  stresses  acting  across  the  elements,  which  is  similar  to  using 
constant  strain  triangular  elements  in  the  FEM.  From  the  FLAG  computed 
stresses,  the  resultant  forces  and  the  points  of  applications  were  computed  for 
both  the  stem  and  heel  sections.  The  resultant  force  and  its  point  of  application 
on  the  stem  are  needed  for  the  structural  design  of  the  stem,  while  the  resultant 
force  and  its  point  of  application  on  the  heel  section  are  required  for  the  global 
stability  of  the  structural  wedge. 

Two  approaches  were  used  to  determine  the  resultant  forces  acting  on  the 
stem  and  heel  sections  from  the  FLAG  computed  stresses.  In  the  first  approach, 
constant  stress  distributions  across  the  elements  were  assumed,  while  in  the 
second  approach,  linearly  varying  stress  distributions  were  assumed.  The  details 
of  the  approaches  are  discussed  in  the  following  subsections. 
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4.1.1  Determination  offerees  assuming  constant-stress  distribution 

The  first  approach  used  to  determine  the  forces  acting  on  the  stem  and  heel 
section  assumed  constant  stress  distributions  across  the  elements,  as  illustrated  in 
Figure  4-1  for  three  of  the  beam  elements  used  to  model  a  portion  of  the  stem. 
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Figure  4-1.  Assumed  constant  stress  distribution  across  elements,  at  time  tj,  used 
to  compute  the  forces  acting  on  the  stem  and  heel  section  in  the  first 
approach 


For  the  assumed  constant  stress  distributions,  the  forces  acting  on  the  top  and 
bottom  nodes  of  each  beam  element,  shown  as  red  dots  in  Figure  4-1,  were 
computed  using  the  following  expressions: 


F‘7=--h 

>,j  2  ' 


top 

ij 


(4-la) 


j^bottom 

iJ 


2 

2 


h  -(T 

I  i,y 


bottom 


(4-lb) 
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where 


FIJ  =  force  acting  on  the  top  node  of  element  i  and  at  time  increment  j 
hi  =  length  of  element  i 

af'j  =  lateral  stress  acting  on  the  top  of  element  i  and  at  time  increment  j 

(7^  J  =  average  lateral  stress  acting  on  element  /  and  at  time  increment  j 

pbouom  _  acting  on  the  bottom  node  of  element  /  and  at  time 
increment  j 

_  lateral  stress  acting  on  the  bottom  of  element  /  and  at  time 
increment  j 

The  total  force  acting  on  the  stem  or  heel  section  Pj  at  time  increment  j  was 
determined  by: 


(4-2) 


4.1.2  Determination  of  forces  assuming  linearly  varying  stress 
distribution 

In  the  second  approach  used  to  determine  the  forces  acting  on  the  stem  and 
heel  sections,  linearly  varying  stress  distributions  across  the  elements  were 
assumed,  as  illustrated  in  Figure  4-2.  This  approach  is  analogous  to  using  a  first- 
order  shape  function  in  the  FEM.  To  apply  this  approach,  the  only  actual  nodes 
considered  were  the  nodes  at  the  top  and  bottom  of  the  wall,  shown  as  red  dots  in 
Figure  4-2.  In  addition  to  these  two  actual  nodes,  the  centers  of  the  elements 
were  treated  as  nodes,  shown  as  green  dots  in  Figure  4-2.  The  magnitudes  of  the 
stresses  on  the  green  nodes  were  assumed  equal  to  the  average  stresses  computed 
by  FLAG  for  the  corresponding  elements.  The  magnitudes  of  the  stresses  acting 
on  the  red  nodes  were  determined  by  linearly  extrapolating  the  stresses  from  the 
neighboring  green  nodes.  In  Figure  4-2,  newly  defined  “elements”  and  “nodes” 
are  numbered  with  primes  (e.g.,  3')  to  distinguish  them  from  the  actual  nodes  and 
elements  shown  in  Figure  4-1.  Additionally,  for  comparison,  the  constant  stress 
distribution  assumed  in  the  previous  approach  is  shown  as  a  dashed  gray  line  in 
Figure  4-2,  as  well  as  the  linearly  varying  stress  distributions  (blue  lines), 
assumed  in  the  current  approach. 
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Figure  4-2.  Assumed  linearly  varying  stress  distribution  across  elements,  at 
time  tj,  used  to  compute  forces  on  stem  and  heel  sections  in  the 
second  approach 


For  the  assumed  linearly  varying  stress  distribution,  the  forces  acting  on  the 
top  and  bottom  nodes  of  the  newly  designated  elements  shown  in  Figure  4-2 
were  computed  using  the  following  expressions. 


(4-3a) 


(4-3b) 


The  total  force  acting  on  a  node  is  the  sum  of  bottom  force  from  the  element 
above  and  the  top  force  from  the  element  below,  i.e.,  and  j ,  respec¬ 

tively.  The  total  force  acting  on  the  stem  and  heel  section  were  computed  using 
the  following  expressions: 


l’,='L(FrT'* 


(4-4) 
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The  resultant  forces  computed  assuming  the  constant  and  linearly  varying 
stress  distributions  were  essential  identical. 


4.1.3  Incremental  dynamic  forces 

In  addition  to  computing  the  total  resultant  forces  acting  on  the  stem  and  heel 
sections,  the  incremental  dynamic  forces  at  time  increment  j  were  computed. 
APy  is  the  difference  between  the  total  resultant  force  Pj  at  time  increment  j  minus 
the  total  resultant  force  prior  to  shaking  (i.e.,  Pj  at  j  =  0,  designated  as  Pstau^'^ 

l^j=Pj-Ps.aac  (4-5) 

Because  Pj  values  computed  assuming  constant  and  linearly  varying  stress 
distributions  were  essentially  identical,  APy  was  computed  only  using  Py  for  the 
constant  stress  distribution. 


4.1.4  Reaction  height  of  forces 


The  points  of  application  of  the  total  and  incremental  dynamic  resultant 
forces  were  computed  for  the  stem  and  heel  sections  in  terms  of  their  vertical 
distances  above  the  base  of  the  wall.  For  the  total  resultant  forces,  the  vertical 
distances  Y  were  computed  using  the  following  relation: 


-y* 


where 


(4-6) 


Yj  -  vertical  distance  from  the  base  of  the  retaining  wall  to  the  point 
of  application  of  the  total  resultant  force  acting  on  the  stem  or  heel 
section  at  time  increment  j 

yi  =  vertical  distance  from  the  base  of  the  retaining  wall  to  the  center 
of  element  / 


The  vertical  distances  AY  from  the  base  of  the  retaining  wall  to  the  points  of 
application  of  APj  acting  on  the  stem  or  heel  section  were  computed  using  the 
following  relation: 


P  Y  -P  Y 

_  J  J  static  static 


(4-7) 
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In  this  equation,  Ysmic  is  the  vertical  distance  from  the  base  of  the  retaining 
wall  to  the  point  of  application  of  the  total  resultant  force  acting  on  the  stem  or 
heel  section  prior  to  the  shaking  (i.e.,  Yj  at  j  =  0). 


4.2  Presentation  and  Discussion  of  Reduced 
Data 

Using  the  procedures  described  in  the  preceding  section,  the  total  and 
dynamic  incremental  resultant  forces  acting  on  the  stem  and  heel  sections  were 
determined  from  the  FLAG  computed  stresses,  as  well  as  the  corresponding 
vertical  distance  above  the  base  at  which  the  resultant  forces  act.  Additionally, 
the  permanent  relative  displacements  of  the  wall  computed  in  the  FLAG  analysis 
are  compared  with  those  predicted  using  a  Newmark  sliding  block-type  analysis 
procedure,  e.g.,  GWROTATE  (Ebeling  and  White,  in  preparation).  Finally,  a 
brief  discussion  is  given  concerning  the  deformed  shape  of  the  wall-soil  system 
at  the  end  of  shaking. 


4.2.1  Total  resultant  forces  and  points  of  action 

The  horizontal  acceleration  a*  and  the  corresponding  dimensionless 
horizontal  inertial  coefficient  kh  at  approximately  the  middle  of  the  backfill 
portion  of  the  structural  wedge  were  computed  during  the  FLAG  analyses,  as 
shown  in  Figure  4-3.  Appendix  B  gives  the  appropriate  sign  convention  related 
to  Oh  and  kh.  In  this  figure,  the  potential  active  and  passive  failure  planes  are 
shown  for  illustration  only.  The  kh  time-histoiy  shown  in  this  figure  is  that  to 
which  reference  is  made  during  the  remainder  of  this  chapter,  unless  otherwise 
noted. 

The  kh  time-history  shown  in  Figure  4-3  is  repeated  in  Figure  4-4a  and  in 
Figure  4-4b  for  reference  purposes  (i.e.,  the  time-histories  in  Figures  4-3, 4-4a, 
and  4-4b  are  identical).  Assuming  constant  stress  distributions  across  the 
elements.  Equation  4-2  was  used  to  compute  the  time-histories  of  the  resultant 
forces  acting  on  the  stem  and  heel  sections  (Pstem  and  Pheei,  respectively), 
presented  as  Figures  4-4c  and  d,  respectively.  Equation  4-6  was  used  to  compute 
time-histories  of  the  vertical  distances  from  the  base  of  the  wall  to  the  point  of 
application  of  the  resultant  forces  acting  on  the  stem  and  heel  sections  (Ys,em  and 
Yheei,  respectively).  These  time-histories  are  presented,  normalized  by  the  height 
of  the  wall  H,  in  Figures  4-4e  and  f.  Finally,  the  time-histories  of  the  resultant 
forces  multiplied  by  the  corresponding  vertical  distances  above  the  base  at  which 
they  act  for  the  stem  and  heel  sections  ((Y-P)sum  and  {Y-P)heei,  respectively)  are 
presented  in  Figures  4-4g  and  h,  respectively. 

Several  interesting  trends  may  be  observed  from  the  time-histories  presented 
in  Figure  4-4.  First,  Ps,em  and  Pheei  are  essentially  at  active  earth  pressure 
conditions  Ka  prior  to  shaking,  which  is  somewhat  unexpected.  The  factor  of 
safety  against  sliding  in  the  static  design  of  the  wall,  per  Gorps  design  procedures 
(refer  to  Appendix  A),  precludes  the  required  movements  to  develop  active 
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Figure  4-3.  Horizontal  acceleration  a/,,  and  corresponding  dimensionless  horizontal  inertial 
coefficient  kfj,  of  a  point  in  the  backfill  portion  of  the  structural  wedge 

conditions  under  static  loading.  Subsequent  FLAG  analyses  of  the  "numerical 
construction”  phase  of  the  wall  will  be  performed  to  more  fully  understand  this 
stress  condition.  Another  observation  concerning  the  lateral  earth  pressures  is 
that  at  the  end  of  shaking,  the  residual  earth  pressures  are  approximately  equal  to 
at-rest  Ko  conditions  for  both  the  stem  and  heel  sections.  Similar  increases  in  the 
earth  pressures  were  found  in  other  studies,  both  numerical  and  laboratory  (i.e., 
centrifuge  and  shake  table),  as  outlined  in  Whitman  (1990).  Additionally,  the 
maximum  value  of  Pheei  is  larger  than  the  maximum  value  of  Pstem,  while  Pstem 
shows  a  much  larger  cyclic  fluctuation  and  at  several  points  dips  slightly  below 
initial  Ka  conditions.  Although  difficult  to  see  from  Figure  4-4  due  to  the  scales, 
Pstem  is  in  phase  with  kh,  while  Pheei  is  out  of  phase  with  kh  (i.e.,  the  peaks  in  one 
time-histoiy  coincide  with  the  troughs  of  the  other  time-history).  Trends  in  the 
relative  magnitude  of  the  resultant  forces  acting  on  the  stem  and  heel  sections, 
and  their  phasing  with  kh,  can  be  more  easily  observed  when  presented  in  terms 
of  the  lateral  earth  pressure  coefficients  K,  The  following  expression  relates  P 
and  A'  at  a  given  time  increment  j: 


2*P 


(4-8) 


where  Kj  is  the  vertical  inertial  coefficient  at  time  increment  j  (assumed  to  be 
zero). 
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Figure  4-4.  Time-histories  of  P,  Y/H  and  YP  for  the  stem  and  heel  sections  (To  convert  kip-feet  to 
Newton-meters,  multiply  by  1,355.8;  to  convert  kips  to  newtons,  multiply  by  4,448) 


Using  this  expression,  K  values  were  computed  for  the  stem  and  heel  sections 
at  the  peaks  and  troughs  during  the  strong  motion  portion  of  the  k/,  time-history 
(i.e.,  5-10  sec,  approximately).  The  K  values  thus  computed  are  plotted  as  func¬ 
tions  of  their  corresponding  absolute  values  of  kh  in  Figure  4-5.  Additionally,  the 
active  and  passive  dynamic  earth  pressure  coefficients  (Kae  and  Kpe,  respec¬ 
tively)  computed  using  the  Mononobe-Okabe  expressions  given  in  Appendix  B 
are  shown  in  Figure  4-5. 
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a.  Entire  range  of  interest 

Figure  4-5.  Comparison  of  lateral  earth  pressure  coefficients  computed  using  the  Mononobe- 
Okabe  active  and  passive  expressions  (yielding  backfill),  Wood  expression  (non¬ 
yielding  backfill),  and  FLAG  (Continued) 

A  portion  of  the  plot  in  Figure  4-5a  is  enlarged  in  Figure  4-5b.  Additionally 
the  kh  time-history  is  given  with  the  peaks  and  troughs  identified  that  correspond 
to  the  computed  K  values.  Several  distinct  trends  may  be  observed  from 
Figure  4-5: 

a.  Kheei  >  K.stem  when  kh<0  (i.e.,  when  ki,  is  directed  toward  the  backfill). 

b.  Ks,em  >  Kheet  whcn  kh>  0  (i.e.,  when  kh  is  directed  away  from  the  backfill). 

c.  The  largest  occurs  when  A:*  >  0  (i.e.,  when  kh  is  directed  away  from 
the  backfill). 
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FLAG  Results 

X  Lateral  E.P.  Coef.  for  heel  section,  kh  away  from  backfill 
O  Lateral  E.P.  Coef.  for  heel  section,  kh  towards  backfill 
X  Lateral  E.P.  Coef.  for  stem,  kh  away  from  backfill 
O  Lateral  E.P.  Coef  for  stem,  kh  towards  backfill 


away  from 
backfill 


towards 

backfill 


b.  Enlargement  of  range  of  FLAG  computed  values 
Figure  4-5.  (Concluded) 

d  The  largest  Kheei  occurs  when  kh<0  (i.e.,  when  kh  is  directed  toward  the 
backfill). 

e.  The  computed  K  values  show  a  general  scatter  around  the  curve  for  the 
Mononobe-Okabe  dynamic  active  earth  pressure  curve. 

The  shape  of  the  Mononobe-Okabe  active  and  passive  dynamic  earth 
pressure  curves  warrant  discussion.  As  kh  increases,  Kae  increases,  while  Kpe 
decreases.  For  the  conditions  examined  (i.e.,  horizontal  backfill,  vertical  wall, 
zero  interface  friction  between  the  structural  and  driving  wedges,  k^  =  0),  Kae  and 
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KpE  reach  the  same  limiting  value.  The  limiting  K  value  occurs  when  the  angles 
of  the  active  and  passive  failure  planes  (which  are  assumed  to  be  planar  in  the 
Mononobe-Okabe  formulation)  become  horizontal;  refer  to  Appendix  B  for 
expressions  for  angles  of  the  failure  planes. 

For  comparison  purposes,  the  earth  pressure  coefficient  for  nonyielding 
backfills  is  also  plotted  in  Figure  4-5.  A  wall  retaining  a  nonyielding  backfill 
does  not  develop  the  limiting  dynamic  active  or  passive  earth  pressures  because 
sufficient  wall  movements  do  not  occur  to  mobilize  the  full  shear  strength  of  the 
backfill,  such  as  is  the  case  with  massive  concrete  gravity  retaining  walls 
founded  on  firm  rock.  Wood  (1973)  developed  a  procedure,  which  was 
simplified  in  Ebeling  and  Morrison  (1992),  Section  5.2,  to  determine  the  lateral 
dynamic  earth  pressures  on  structures  with  nonyielding  backfills.  The  following 
expression  is  from  Ebeling  and  Morrison  (1992),  Equation  68: 

Fsr=yH^-h  (4-9) 


where 


Fsr  =  lateral  seismic  force  component 
y  =  unit  weight  of  the  soil 

By  treating  Fsr  as  the  dynamic  incremental  force,  the  equivalent  earth 
pressure  coefficient  was  computed  by  substituting  Fsr  into  Equation  4-8  for  P  and 
adding  Ko  to  the  result.  The  resulting  curve,  shown  in  Figure  4-5,  will  likely  be  a 
conservative  upper  bound  of  the  earth  pressures  that  will  occur  on  the  heel 
section  of  a  cantilever  wall.  However,  a  more  probable  upper  bound  is  that 
formed  by  a  line  drawn  from  Ko  pressure  for  kh  =  0  and  the  intersection  of  Kae 
and  KpE  at  their  limiting  values.  Further  FLAG  analyses  will  be  performed  to 
verify  this  hypothesized  upper  bound. 

Similar  to  the  trends  in  Pstem  and  Pheeh  YIH  for  the  stem  and  heel  sections  (i.e., 
YIHheei  and  YIHstem,  respectively)  also  show  increasing  trends  as  the  shaking 
progresses,  with  YIHstem  having  greater  cyclic  fluctuation  than  YIHheeh  Of 
particular  note  is  that  YIHstem  is  out  of  phase  with  both  kh  and  Pstem,  while  YIHheei 
is  out  of  phase  with  kh,  but  in  phase  with  Pheeh  As  a  result  of  the  phasing,  {YP)heei 
has  considerably  larger  cyclic  fluctuations  and  peak  value  than  {YP)stem- 

The  magnitudes  of  Ystem  and  Yheei  are  directly  related  to  the  distribution  of 
stresses  along  the  stem  and  heel  sections,  respectively.  The  stress  distributions, 
resultant  forces,  and  deformed  shape  of  the  cantilever  wall  corresponding  to 
maximum  values  of  Pstem,  Pheei,  {YP)stem,  and  (YP)heei  are  shown  in  Figure  4-6, 
where  the  maximum  values  for  Pheei,  (YP)stem,  and  (YP)heei  all  occur  at  the  same 
instant  in  time.  The  maximum  value  of  Pstem  occurs  while  kh>0  (i.e.,  kh  is 
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directed  away  from  the  backfill);  refer  to  the  kh  time-histoiy  in  Figure  4-5b.  The 
relatively  triangular-shaped  stress  distributions  on  the  stem  and  heel  sections 
shown  in  Figure  4-6a  are  characteristic  of  those  occurring  at  the  peaks  in  the  kh 
time-histoiy.  The  points  of  action  of  the  resultant  forces  on  stem  and  heel 
sections  are  approximately  equal  to  those  prior  to  the  start  of  the  shaking. 

The  maximum  values  of  Pheeh  {YP)stem,  and  {YP)heei  occur  while  kh<0  (i.e., 
kh  is  directed  toward  the  backfill);  again,  refer  to  the  kh  time-history  in  Fig¬ 
ure  4-5b.  The  relatively  uniform  stress  distributions  on  the  stem  and  heel 
sections  shown  in  Figure  4-6b  are  characteristic  of  those  occurring  at  the  troughs 
in  the  kh  time-history.  The  points  of  action  of  the  resultant  forces  on  stem  and 
heel  sections  are  approximately  at  midheight  of  the  wall,  and  therefore  higher 
than  the  static  values. 


4.2.2  Ratio  of  total  resultant  forces  and  points  of  action 

Given  the  difference  in  the  phasing  of  the  stem  and  heel  section  time- 
histories,  it  is  only  of  limited  value  to  examine  the  ratios  of  the  respective  time- 
histories  for  the  stem  and  heel  sections,  which  are  presented  in  Figure  4-7.  For 
example,  from  examination  of  Figure  4-7b,  PsteJ Pheeh  one  might  conclude  that 
the  maximum  value  for  Pstem  is  larger  than  the  maximum  value  for  Pheeh  which  is 
clearly  not  the  case,  as  may  be  seen  from  Figure  4-4c  and  d.  Similar  misinter¬ 
pretations  can  be  made  in  relation  to  the  ratios  of  the  other  time-histories 
presented  in  Figure  4-7. 


4.2.3  Incremental  resultant  forces  and  points  of  action 

As  an  alternate  to  presenting  the  total  resultant  force  of  the  lateral  earth 
pressures.  Seed  and  Whitman  (1970)  expressed  the  resolved  lateral  earth 
pressures  in  terms  of  a  static  active  resultant  {P static)  and  a  dynamic  incremental 
resultant  (AP),  as  discussed  previously  in  Section  4.1.3.  Seed  and  Whitman's 
(1970)  procedure  is  outlined  in  Ebeling  and  Morrison  (1992),  Section  4.2.2. 

Using  Equations  4-5  and  4-7,  time-histories  for  AP  acting  on  the  stem  and 
heel  sections  and  their  corresponding  height  AT  above  the  base  of  the  wall  were 
determined  from  the  FLAG  computed  stresses.  The  computed  time-histories  for 
AP  and  AF-AP  are  shown  in  Figure  4-8.  Pstatic  is  included  in  Figure  4-8c  and  d 
for  reference  purposes  only.  The  computed  time-histories  of  AF  were  not 
included  in  this  figure  because  the  extreme  ranges  of  the  AF  values  and  the 
erratic  characteristics  of  the  time-histories  made  them  impossible  to  present  in  an 
intelligible  manner. 

In  addition  to  the  results  computed  from  the  FLAG  stresses,  the  simplified 
procedure  proposed  by  Seed  and  Whitman  (1970),  and  outlined  in  Ebeling  and 
Morrison  (1992),  Section  4.2.2,  was  used  to  compute  AP,  AF,  and  AP-AF,  i.e.,: 
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Figure  4-8.  Time-histories  of  AP  and  AY  AP  for  the  stem  and  heel  sections  (To  convert  kips  to 
newtons,  multiply  by  4,448;  to  convert  kip-feet  to  newton-meters,  multiply  by 
1,355.8) 


The  validity  of  these  expressions  is  limited  to  specific  values  of  (]),  0,  6,  and  P 
(the  refer  to  Appendix  B  for  definitions  of  these  variables).  More  general 
expressions  for  AP  and  AY,  which  are  likely  more  applicable  to  most  Corps 
projects,  are  given  in  Ebeling  and  Morrison  (1992),  Section  4.2.2. 

Setting  kh  equal  to  the  maximum  peak  value  (i.e.,  =  0.32),  seed  and  whumam 

tAY Seed  and  Whitman,  and  {AY •b^)seed  and  Whitman  Were  computed  usittg  Equations  4-1  Oa, 
b,  and  c;  and  IaP seed  and  whitman  and  {AY  AP)seed  and  whitman  are  presented  in  Figure  4-8. 
As  may  be  observed  in  Figure  4-8c,  /AP seed  and  whitman  is  very  close  to  the  maximum 
value  of  M^stem  from  the  FLAG  results.  However,  this  may  be  coincidental  given 
that  /APstem  does  not  coincide  with  the  maximum  peak  in  the  time-history. 

/aP Seed  and  Whitman  IS  Icss  than  maximum  value  of  lAPheeh  which  is  associated  with  a 
high  frequency  spike.  Most  notable  of  the  observed  trends  in  Figure  4-8  is  that 
the  maximum  value  of  (A7-AP),,,.;„  is  considerably  less  than  the  maximum  value 
of  (AT'AP)/,^^,/  and  {/AY'/AP^seedandWhuman^  Ih  regard  to  the  comparison  of  the 
maximum  values  of  (AT'AP)^,^;,,  and  {AY-AF)heeU  the  difference  is  due  largely  to 
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the  phasing  of  the  respective  AF  and  A/*  time-histories,  where  and  APs,em 

are  out  of  phase  with  each  other,  while  AYheei  and  APneei  are  in  phase. 

Analogous  to  Figure  4-6,  Psiatic  and  AP  for  the  stem  and  heel  sections  are 
shown  in  Figure  4-9  superimposed  on  the  stress  distributions  and  deformed  shape 
of  the  cantilever  wall  at  times  corresponding  to  maximum  values  of  Psiem,  Pheei, 
{YP)s,em,  and  {YP)heei-  Of  particular  note  are  heights  above  the  base  at  which  AP 
on  the  stem  and  heel  sections  act  (i.e.,  AYsum  and  AYf^ei,  respectively).  For  the 
condition  where  Pstem  is  maximum,  AY^tem  is  approximately  equal  to  Ystaac,  which 
is  about  0.3  Pf.  However,  AYhed  is  closer  to  0.5-H.  For  the  condition  when  Phed, 
(YP)s,em,  and  (YP)heei  are  maximum,  AYs^m  is  just  over  0.75  H,  while  AY/^ei 
remains  at  approximately  0.5  //.  Figure  4-9  also  illustrates  the  phasing  of  the 
various  parameters.  For  the  condition  where  Ps^m  is  maximum,  APsiem  is  large, 
while  AYstem  and  APw  are  relatively  small.  For  the  condition  where  Pheei, 
(YP^iem,  and  (YP)heei  are  maximum,  APsiem  is  small,  while  AYs,em  and  APheei  are 
relatively  large. 


4.2.4  Permanent  relative  displacement  of  the  wall 

Using  an  acceleration  time-history  computed  by  FLAG  at  approximately 
middepth  of  the  backfill  and  located  near  the  free-field  boundary  in  the  FLAG 
model,  a  Newmark  sliding  block-type  analysis  was  performed  on  the  structural 
wedge.  This  emalysis  was  similar  to  those  performed  using  GWROTATE 
(Ebeling  and  White,  in  preparation)  with  the  SHAKE  computed  time-histories, 
which  are  presented  in  Appendix  G.  The  results  from  the  sliding  block-type 
analysis  were  compared  to  wall  movements  computed  by  FLAG. 

The  FLAG  computed  free-field  acceleration  time-history  used  in  the  sliding 
block  analysis  is  shown  in  Figure  4-lOa.  Also  shown  in  Figure  4-lOa  is  the 
acceleration  time-history  of  the  structural  wedge,  as  determined  by  sliding  block 
einalysis.  Relative  movement  of  the  soil  and  the  structural  wedge  occurs  when 
the  free-field  acceleration  exceeds  the  maximum  transmissible  acceleration 
(N*-g)  of  the  structural  wedge,  where  N*-g  was  determined  to  be  0.22g  in 
Appendix  G.  Figure  4-1  la  shows  an  enlargement  of  the  region  where  the  two 
acceleration  time-histories  differ;  Figure  4-1  lb,  c,  and  d  show  the  progression  of 
the  steps  used  to  compute  the  permanent  relative  displacement  using  a  Newmark 
sliding  block-type  analysis.  Figure  4- 10b  shows  a  comparison  of  the  FLAG 
computed  permanent  relative  displacement  time-history  with  that  computed  in 
the  sliding  block  analysis.  The  magnitude  and  the  occurrence  times  of  the  slips 
are  very  similar  for  both  analyses.  However,  the  FLAG  results  have  what  is 
believed  to  be  a  “numerical  creep”  prior  to  approximately  7  sec  and  after 
approximately  20  sec.  This  numerical  creep  is  likely  due  to  precision  error  that 
occurs  when  the  displacement  experienced  during  one  time  increment  is 
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Figure  4-9.  Stress  distributions,  static  and  incremental  dynamic  resultant  forces  on  the  stem  and 
heel  sections  at  times  corresponding  to  the  following;  (a)  maximum  value  for  Pstem.  and 
(b)  the  maximum  values  for  P^eei.  {YP)stem,  and  {y  P)heei  (To  convert  feet  to  meters, 
multiply  by  0.3048;  to  convert  psf  to  pascals,  multiply  by  47.88;  to  convert  kips 
to  newtons,  multiply  by  4,448) 
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Figure  4-10.  Comparison  of  the  permanent  relative  displacements  predicted  by 
a  Newmark  sliding  block-type  analysis  and  by  FLAG:  (a)  k/,  time- 
history  used  in  the  analysis,  (b)  sliding  block  and  FLAG  results 
(To  convert  inches  to  millimeters,  multiply  by  25.4) 


relatively  small  compared  with  the  cumulative  total  displacement.'  In  cases 
where  higher  accuracy  is  required  in  the  cumulative  relative  displacement,  rather 
than  just  the  magnitude  of  the  slip  increments,  the  double  precision  version  of 
FLAG  should  be  used. 


4.2.5  Deformed  grid  of  the  wall-soil  system,  post  shaking 

Figure  4-12  shows  the  deformed  grid  of  the  wall-soil  system  after  the 
completion  of  the  earthquake  shaking,  magnified  by  a  factor  of  ten.  Two 
interesting  observations  may  be  made  concerning  this  figure.  First,  the  interface 
between  the  structural  and  driving  wedges  does  not  appear  to  be  vertical,  as  is 
often  assumed  in  simplified  analysis  procedures.  The  practical  significance  of 
this  needs  to  be  explored  further. 

A  second  observation  is  that  as  opposed  to  one  distinct  failure  plane  running 
from  the  heel  of  the  wall  up  through  the  backfill  to  the  right,  multiple  failure 
planes  or  shear  bands  are  observed.  The  range  of  the  shear  bands  is  highlighted 


*  N.  McCullough,  2002,  Personal  communication  regarding  FLAG  modeling.  Department  of  Civil, 
Construction,  and  Environmental  Engineering,  Oregon  State  University,  Corvallis,  OR. 


Chapter  4  FLAG  Data  Reduction  Discussion  of  Results 


47 


Time  (sec) 


Figure  4-1 1 .  Results  from  the  Newmark  sliding  block-type  analysis  of  the  structural  wedge; 

(a)  acceleration  time-histories  of  the  soil  and  structural  wedge,  (b)  absolute  velocity 
time-histories  of  the  soil  and  structural  wedge,  (c)  relative  velocity  of  the  soil  and 
structural  wedge,  and  d)  permanent  relative  displacement  of  the  structural  wedge 
(To  convert  inches  to  millimeters,  multiply  by  25.4) 
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Figure  4-12.  Deformed  grid  of  the  wall-soil  system,  post  shaking,  magnification 
X  10:  (a)  without  markups,  (b)  with  markups 


in  Figure  4- 12b.  The  formation  of  multiple  failure  planes  is  consistent  with 
Equation  B-5a,  which  predicts  the  angle  of  inclination  of  the  failure  plane  to 
flatten  as  kh  increases.  Accordingly,  it  would  be  expected  that  a  range  of  failure 
to  plane  would  form  during  the  random  shaking  of  earthquake  motions.  Fig¬ 
ure  4-13  presents  model  wall  tests  performed  on  a  shake  table.  Multiple  failure 
planes  are  clearly  identifiable  in  these  photographs.  However,  the  formation  of 
multiple  failure  planes  is  not  always  observed  in  laboratory  studies.  For  dilative 
soils,  the  formation  of  one  failure  plane  may  be  favored  over  the  formation  of 
multiple  failure  planes. 


4.3  Conclusions 

In  the  preceding  sections,  the  procedures  used  to  reduce  the  FLAG  data  were 
presented  for  this,  the  first  FLAG  analysis  of  a  Gorps  type  retaining  wall. 
Additionally,  the  reduced  data  were  presented  and  preliminary  interpretation 
given.  Further  studies  are  required  to  determine  whether  the  trends  identified  in 
the  first  phase  of  this  study  are  general  in  nature,  or  limited  to  certain  wall 
configurations,  assumed  wall  and  soil  properties,  and  soil  conditions. 
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Figure  4-1 3.  Shake  table  tests  performed  on  scale  models  of  retaining  wall.  Note 
the  multiple  failure  planes  in  the  backfill  (from  Aitken,  Elms,  and 
Berrill  1982) 
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Appendix  A 

Static  Design  of  the  Cantilever 
Retaining  Wall 


A.1  Introduction 

This  appendix  illustrates  the  sizing  and  structural  design  for  usual  loading 
(i.e.,  static)  of  the  20-ft-high*  cantilever  retaining  wall  that  is  analyzed  dynamic¬ 
ally  in  the  main  body  of  this  report.  The  wall  design  is  performed  in  two  stages. 
The  first  stage  consists  of  sizing  the  wall  to  satisfy  global  stability  requirements, 
in  general  accordance  with  Engineer  Manual  (EM)  1 1 10-2-2502  (Headquarters, 
U.S.  Army  Corps  of  Engineers  (HQUSACE)  1989).^  The  global  stability 
requirements  are  expressed  in  terms  of  the  factor  of  safety  against  sliding;  the 
factor  of  safety  against  bearing  capacity  failure;  and  the  percentage  of  the  base 
area  in  compression,  with  the  latter  quantifying  the  overturning  stability  of  the 
wall. 

The  second  stage  of  the  wall  design  entails  the  dimensioning  of  the  com¬ 
ponents  of  the  concrete  wall  (i.e.,  stem  and  base  slab,  toe  and  heel  elements)  and 
detailing  of  the  reinforcing  steel,  in  general  accordance  with  EM  1 1 10-2-2104 
(HQUSACE  1992).  Each  of  the  three  structural  elements  is  designed  by  the 
strength-design  method  as  a  cantilever,  one-way  slab  for  flexure  and  shear 
loadings.  The  design  loads  are  those  determined  in  the  first  design  stage,  with 
appropriate  dead,  live,  and  hydraulic  load  factors  applied,  thus  ensuring  that  the 
Corps  strength  and  serviceability  requirements  are  satisfied. 

The  notation  used  in  this  appendix  is  generally  consistent  with  that  used  by 
the  Corps.  However,  both  structural  and  geotechnical  calculations  are  presented 
in  this  appendix,  and  the  notation  between  the  two  disciplines  is  somewhat 
ambiguous  (e.g.,  (|)  is  used  in  both  concrete  design  and  geotechnical  engineering, 
with  very  distinct  meanings).  Accordingly,  rather  than  giving  a  compiled  list  of 
veuiable  definitions  used  in  this  appendix,  variables  are  defined  mathematically 
and  illustratively  as  they  are  presented  in  the  calculations. 


*  Conversion  factors  for  non-SI  units  of  measurement  used  in  this  appendix  are  given  in 
Table  A-1. 

^  References  cited  in  this  appendix  are  included  in  the  References  section  at  the  end  of 
the  main  text. 
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A.2  Stage  1 :  Sizing  of  the  Cantiiever  Retaining 
Wail 

As  stated  previously,  the  first  design  stage  consists  of  sizing  the  cantilever 
wall  such  that  global  stability  requirements  are  satisfied  (i.e.,  sliding,  overturn¬ 
ing,  and  bearing  capacity),  in  general  accordance  with  EM  1 1 10-2-2502 
(HQUSACE  1989).  The  structural  wedge  of  the  proposed  wall  and  backfill  is 
shown  in  Figure  A-1,  as  well  as  the  backfill  and  foundation  material  properties. 

To  assess  the  global  stability  of  the  wall,  the  external  forces  and  correspond¬ 
ing  points  of  action  acting  on  the  structural  wedge  need  to  be  determined.  The 
external  forces  include  the  resultant  of  the  lateral  earth  pressure  and  reactionary 
forces  acting  along  the  base  of  the  wall.  However,  before  the  reactionary  forces 
can  be  determined,  the  weights  and  centers  of  gravity  of  the  concrete  and  soil 
composing  the  structural  wedge  are  required. 


Backfill;  medium-dense  cohesionless  compacted  fill,  Yn,  =  125  pcf,  (j>  =  35° 
Foundation:  thin  layer  of  compacted  fill  overlying  natural  deposit  of  dense 
cohesionless  soil,  Yn,  =  125  pcf,  (j)  =  40° 

Reinforced  concrete:  y  =  1 50  pcf,  fe  =  4  ksi,  fy  =  48  ksi 
Hydraulic  factor:  1 .3 


Figure  A-1 .  Structural  wedge  of  proposed  wall,  and  backfill  and  foundation  properties 


A2 


Appendix  A  Static  Design  of  the  Cantilever  Retaining  Wall 


To  simplify  the  determination  of  the  weight  and  center  of  gravity  of  the 
structural  wedge,  it  is  divided  into  subsections  having  uniform  unit  weights  and 
simple  geometries,  as  shown  in  Figure  A-2.  Once  the  weights  and  centers  of 
gravity  of  each  subsection  are  determined,  the  weight  and  center  of  gravity  of  the 
entire  structural  wedge  are  easily  determined,  as  illustrated  in  the  accompanying 
sketches. 


Backfill: 


forces  and  centroids  of  structural 
wedge 


Ws  =  125pcfx8'x  18’ 

=  18,000  lbs 

Axtoe-Ws  =  18,000  lbs  X  9' 
=  162,000  ft-lbs 


00 


ToeH" 


Axtoe  =  9  ft 
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Wc  stem  2  =  1 50pcf  X  0.5’  X  1 8’  X  0.5  Wc  stem  I  =  1 50pcf  X  1 .5'  X  1 8’ 

=  675  lbs  =  4050  lbs 

Ax,oe-Wc  stem  2  =  675  Ibs  X  3.33’  Ax,oe- Wc  stem  I  =  4050  Ibs  X  4.25’ 

=  2250ft-lbs  =  17,212.5  ft-lbs 


Wcbase=150pcfx2’x  13’ 

=  3900  lbs 

Ax,oe-Wcbase  =  3900  Ibs  X  6.5’ 
=  25,350  ft-lbs 
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Weight  of  entire  wall; 


^Wc,i  =  3900  lbs  +  4050  lbs  +  675  lbs 
i=l 

=  8625  lbs 


Centroid  of  entire  wall: 

3 

^Wc.i  •  Ax,oe,i  =  25,350  ft-lbs  +  17,212.5  ft-lbs  +  2250  ft-lbs 
i=  1 

=  44,812.5  ft-lbs 
3 

Wj^  i  •  AXtoe,  i 

- =  5.2  ft  from  the  toe 

Zw., 

i=  1 


In  a  CTWALL  analysis  the  retaining  wall  system  is  divided  into  two  or  three 
wedges:  the  structural  wedge;  the  driving  wedge;  and  the  resisting  wedge,  when 
present  (Pace  1994).  Figure  A-1  shows  the  structural  wedge,  which  is  defined  by 
the  outline  of  the  cantilever  retaining  wall.  The  lateral  extent  of  the  structural 
wedge  is  defined  by  imaginary  vertical  sections  made  through  the  heel  of  the 
wall  and  the  toe  of  the  wall.  The  soil  mass  contained  within  this  region  is  also 
considered  part  of  the  structural  wedge.  The  driving  soil  wedge  on  the  retained 
soil  side  (to  the  right  of  the  Figure  A-1  structural  wedge  and  not  shown  in  this 
figure)  generates  earth  pressure  forces  tending  to  destabilize  the  structural  wedge. 
No  resisting  wedge  is  present  to  the  left  of  the  Figure  A-1  structural  wedge  in  this 
case. 

The  general  wedge  method  of  analysis  (EM  1 1 10-2-2502  (HQUSACE 
1989))  is  used  to  calculate  the  lateral  earth  pressure  force  acting  on  the  structural 
wedge.  CTWALL  (Pace  1994)  performs  the  sliding  and  overturning  stability 
analyses  during  execution.  Forces  computed  during  the  CTWALL  analysis  of  the 
user-specified  wall  geometry  may  be  used  to  compute  the  factor  of  safety  against 
a  foundation  bearing  failure.  This  computation  is  made  external  to  the  CTWALL 
execution. 

A  sliding  stability  analysis  using  CTWALL  (Pace  1994)  follows  Corps 
design  criteria  (EM  1 1 10-2-2502).  A  sliding  stability  analysis  is  conducted  of  the 
driving  wedge,  structural  wedge,  and  resisting  wedge  (when  present)  to 
determine  a  common  factor  of  safety  against  sliding  for  the  entire  retaining 
structural  system  (of  the  three  wedges).  The  procedure  is  a  limit  equilibrium 
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procedure  and  is  iterative  in  nature.  Critical  failure  angles  for  the  driving  side 
and  resisting  side  (when  present)  potential  planar  slip  planes  are  sought  for  an 
assumed  sliding  factor  of  safety  value.  The  resulting  earth  forces  are  summed.  If 
the  sum  is  zero,  the  system  is  in  equilibrium  and  the  critical  sliding  factor  of 
safety  value  has  been  found. 

An  overturning  stability  analysis  using  CTWALL  (Pace  1994)  follows  Corps 
design  criteria  (EM  1 110-2-2502).  The  program  will  calculate  a  percentage  of  the 
base  in  compression  and  report  the  location  of  the  resultant  of  all  forces  applied 
to  the  structural  wedge.  The  program  will  iterate,  if  necessary,  to  find  the  percent 
of  the  base  in  compression.  Tension  along  the  base  is  not  allowed  for  these  soil- 
founded  walls.  A  linear  base  pressure  distribution  assumption  is  made  by 
CTWALL  (for  the  portion  of  the  base  in  compression)  in  conjunction  with  the 
limit  equilibrium  procedure  of  analysis.  The  overturning  analysis  checks  that  a 
user-specified  wall  configuration  satisfies  Corps  design  criteria.  These  criteria 
are  expressed  in  terms  of  base  area  in  compression  (versus  computing  the  ratio  of 
the  stabilizing  moment  about  the  toe  of  the  wall  divided  by  the  overturning 
moment).  The  overturning  analysis  uses  an  approximation  for  at-rest  earth 
pressure  forces  in  the  driving  wedge  analysis  to  determine  the  earth  force  the 
driving  wedge  (i.e.,  the  retained  soil)  exerts  on  the  structural  wedge.  Per  EM 
1 1 102-2502  (section  3-11),  the  shear  mobilization  factor  (SMF)  is  set  equal  to 
2/3.  The  SMF  and  factor  of  safety  are  inverses  of  each  other.  An  SMF  of  2/3  is 
equivalent  to  a  factor  of  safety  equal  to  1.5.  The  determination  of  the  lateral  earth 
forces  and  pressures  are  illustrated  in  the  following  equations  in  simplified  hand 
computations  (and  were  used  to  check  the  CTWALL  analyses  for  the  final  wall 
configuration): 


SMF  =  2/3  =  1/1.5 

tan  ((l)'niob)  =  SMF  •  tan 

=  .  tan  (35°) 

=>  <t>'inob  =  25° 


In  accordance  with  EM  1 1 10-2-2502,  the  angle  of  friction  alongside  the 
structural  wedge  defined  by  the  vertical  section  through  the  heel  of  the  wall  is 
assumed  equal  to  zero  (i.e.,  6  =  0°).  Thus,  for  the  case  being  considered  (i.e., 
homogenous  backfill),  the  wedge  analysis  procedure  reverts  to  the  classical 
Rankine  procedure,  which  is  illustrated  in  the  following  sketches. 
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k.  =  tan=(45»-  — ) 


25° 

=  tan^(45°- - ) 

2 

=  0.41 

^  h  “  Ymoist  ■  H  •  kh 

=  125pcf-20'-0.41 
=  1013.8  psf 


Fh. 


Static 


a’h-H 

2 

^—  •1013.8  psf -20' 
2 

10,137.5  lbs 


20  ft 


Aytoe 


X 


Aytoe  = 


H 

3 

20' 

3 

6.667  ft 


Fh,  static  X  Aytoe  =  10,137.5  lbs  •  6.667' 
=  67,583  ft-lbs 


The  remaining,  yet-to-be-determined,  forces  acting  on  the  structural  wedge 
are  the  reactionary  shear  and  normal  forces  (i.e.,  T  emd  If,  respectively)  acting  on 
the  base  of  the  wall.  Figure  A-3  shows  a  free  body  diagram  of  the  structural 
wedge  with  the  known  and  unknown  forces  identified.  As  illustrated  in  the 
following  equations,  the  magnitude  of  T  and  If  are  determined  by  summing  the 
forces  in  the  horizontal  and  vertical  directions,  respectively,  while  the  point  at 
which  A*  acts  is  determined  by  summing  the  moments  around  the  toe. 


+t  Zf'-" 

=  N'-Wc-Ws 
=>  N'=  Wc  +  Ws 

=  8625  lbs  +  18,000  lbs 
=  26,625  lbs 


=  10,137.5  lbs 
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With  all  external  forces  acting  on  the  structural  wedge  determined,  the  global 
stability  of  the  wall  is  assessed  by  computing  the  factor  of  safety  against  sliding; 
the  percentage  of  the  base  area  in  compression  (overturning  stability);  and  the 
factor  of  safety  against  bearing  capacity  failure. 


A.2.1 .  Factor  of  safety  against  sliding 

EM  1 1 10-2-2502  requires  the  following  minimum  factors  of  safety  against 
sliding  (FSsiiding)  for  cantilever  retaining  walls: 

a.  1.5:  usual  loadings. 

b.  1 .33 :  short-duration  (unusual)  loads,  such  as  those  that  might  occur 
during  high  winds,  construction  activities,  or  the  Operational  Basis 
Earthquake  (OBE). 

c.  1.1:  extreme  loads,  such  as  those  that  might  occur  during  the  Maximum 
Design  Earthquake  (MDE). 

The  FSsiiding  for  usual  loadings  is  computed  eis  follows: 

full  —  N  •  tan(5|)3se)  FSsiiding  ~ 

5base  —  ^l^base 

=  26,625  lbs  tan(35°) 

=  18,643  lbs 


Tuit 

T 

18,643  lbs 
10,137.5  lbs 

1.84  >1.5  okay 


Note:  It  is  assumed  that  the  wall  sits  on  a  thin  layer  of  compacted  fill.  The 
internal  friction  angle  of  this  thin  layer  is  assumed  equal  to  that  of  the  compacted 
backfill  (i.e.,  (j)  =  35°). 


A.2.2  Percentage  of  the  base  area  in  compression  (i.e.,  overturning 
stability) 

The  global  stability  of  cantilever  retaining  walls  to  overturning  is  quantified 
by  the  percentage  of  the  base  area  in  compression.  EM  1 1 10-2-2502  requires  the 
following  minimum  percentages  cantilever  walls  on  soil  foundations: 

a.  100%:  usual  loadings. 

b.  75%:  unusual  loadings. 

c.  Resultant  within  base:  extreme  loadings. 

The  percentage  of  the  base  area  that  is  in  compression  for  usual  loadings  is 
computed  as  follows: 
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Cmin^O  /.  100%  base  area  in  compression:  okay 

Note:  Per  Corps  design  criteria,  a  linear  effective  base  pressure  is  assumed. 

A.2.3  Factor  of  safety  against  bearing  capacity  failure 

EM  1 1 10-2-2502  requires  the  following  minimum  factors  of  safety  against 
bearing  capacity  failure  FSbc  for  cantilever  retaining  walls: 

a.  3.0:  usual  loadings. 

b.  2.0:  unusual  loadings. 

c.  >1 :  extreme  loadings. 

EM  1 1 10-2-2502  provides  the  following  expression  for  the  normal 
component  to  the  base  of  the  structure  of  the  ultimate  bearing  capacity  for  strip 
footings: 

Q  =  B  [(^cd-^ci-^cf^cg-C-Nc)  +  - -] 
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However,  for  the  wall  being  analyzed,  only  the  last  term  is  nonzero,  and  thus, 
this  expression  reduces  to: 


Q=  B 


The  FSbc  for  usual  loadings  is  computed  as  follows: 


B  =  B-2e 
=  13'  -  2  (1.27) 
=  10.46' 

5.=  (i-  4)' 

<l> 

=  (l.i^) 


40° 


=  0.23 


g  -  B  - - - - 


-  tan-‘  r  t 

^  26,625  lbs  ^ 

=  20.8° 

^yd  ~  “  ^yg  “  1 


FSbc- 


Q. 

N' 


=  10.46' 

=  147,612  lbs 


(1  •  0.23  •  1  •  1  •  10.46'  •  125pcf  •  93.69) 
2 


^  147,612  lbs 
26,625  lbs 

=  5.54  okay 


Note:  Although  the  wall  sits  on  a  thin  layer  of  compacted  fill  having  (j)  =  35°, 
this  layer  will  have  little  influence  on  bearing  capacity,  as  the  failure  surfaces 
will  primarily  pass  through  the  underlying  denser  natural  deposit  (<|)  =  40°).  This 
is  contrary  to  the  sliding  mode  of  failure  where  the  thin  layer  of  compacted  fill 
has  significant  influence,  and  accordingly,  (j)  =  35°  was  used  to  compute  FSsiiding- 

A.3  Stage  2:  Structural  Design  of  Concrete 
Cantilever  Retaining  Wall 

As  stated  in  the  introduction  to  this  appendix,  the  second  stage  of  the  wall 
design  entails  the  structural  design  of  the  concrete  wall,  to  include  the  dimension¬ 
ing  of  the  concrete  base  slab  (the  toe  and  heel  elements)  and  stem,  and  the  detail¬ 
ing  of  the  reinforcing  steel.  All  reinforced-concrete  hydraulic  structures  must 
satisfy  both  strength  and  serviceability  requirements.  In  the  strength  design 
method,  this  is  accomplished  by  multiplying  the  service  loads  by  appropriate 
load  factors  and  by  a  hydraulic  factor.  The  hydraulic  factor  is  used  to  improve 
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crack  control  in  hydraulic  structures  by  increasing  reinforcement  requirements, 
thereby  reducing  steel  stresses  at  service  loads.  The  service  loads  are  those 
determined  in  the  first  design  stage  presented  previously. 

Each  of  the  three  structural  elements  is  designed  as  a  cantilever,  one-way 
slab  for  flexure  and  shear  loadings,  in  general  accordance  with  EM  1 1 10-2-2104 
(HQUSACE  1992),  with  the  exception  of  the  grade  steel  of  the  reinforcing  bars. 
The  following  example  uses  Grade  40  steel,  which  is  in  line  with  the  design 
examples  in  EM  1 1 10-2-2502  (HQUSACE  1989).  However,  EM  1 1 10-2-2104 
(Section  2-2)  recommends  that  Grade  60  reinforcing  steel  be  used.  Four-inch 
cover  is  used  in  the  example,  per  EM  1 1 10-2-2104  (Section  2-6).  Where  design 
criteria  is  absent  in  EM  1 1 10-2-2104  and  EM  1 1 10-2-2502,  ACl  3 1 8  (American 
Concrete  Institute  2002)  was  used.  Figure  A-4  shows  the  structural  wedge  and 
the  externally  imposed  stresses  determined  in  the  first  design  stage.  Also  shown 
in  Figure  A-4  are  the  critical  locations  for  evaluating  shear  and  bending  moment 
for  the  stem,  heel,  and  toe  elements. 


k -  13  ft - H 


Figure  A-4.  Critical  locations  for  shear  and  bending  moment,  per  ACl 
318  and  EM  1110-2-2502 
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A.3.1  Moment  capacity  of  the  stem 


The  stem  is  analyzed  as  being  singly  reinforced  with  the  critical  section  for 
moment  capacity  being  at  the  base  of  the  stem,  as  illustrated  below. 


=  95,147  ft-lbs  =  1 142  in-kips  Mstem 


Minimum  Effective  Depth: 


dd  = 


2.4956  •  M„ 


2.4956  •  1 142  in-kips 
IT' 


=  15.4"  <  d  okay 


Minimum  Required  Reinforcing  Steel: 


ku=  1 


M„  +  P„-(d-h/2) 
0.425  ■  fc-b  d^ 


=  1  - 


1142  in-kips 

0.425- 4  ksi- 12" -(19.5")^ 


0.85  •  f  c-ku-b-d 


0.85-4  ksi  0.0765  12"19.5" 
48ksi 

=  1 .27  in^  per  1  ft  width  of  wall 


=  0.0765 


Use  #11  @9"  c-c  {conservative) 

As  =  2.08  in^  per  ft  of  wall 
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Check  Ductility  Requirements: 


Pmin  S  P  S  Pmax 


200 

48,000 


=  0.00417 


Pmax  0.375'Pb 

=  0.375  •  0.0388 
=  0.0146 


^  2.08  in^ 

12"-  19.5" 

=  0.0089 


Pmin  S  P  S  Pmax  okcty 

A.3.2  Shear  capacity  of  the  stem 


The  critical  section  for  shear  in  the  stem  is  taken  as  19.5"  above  the  interface 
of  the  base  and  stem,  where  19.5"  is  d  at  the  base  of  the  stem.  However,  the  d  at 
the  critical  section  is  only  19",  due  to  the  taper  of  the  wall. 


— 1\  19"  [« — 
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(l)-Vs>1.3-(V„-fVe) 

1.3  (Vu  -  (l)  Vc)  =  1.3  (1 1,681  lbs  -  24,514  lbs) 

<0  okay 

A.3.3  Moment  and  shear  capacity  of  the  heel 

The  heel  is  analyzed  as  being  singly  reinforced,  with  the  steel  along  the  top 
face  and  4  in.  coverage.  The  critical  section  for  both  moment  and  shear  capacity 
in  the  heel  is  at  the  interface  of  the  heel  and  the  stem. 


4V^MheeI  =  0 

=>  4'  •  (18,000  lbs  +  2400  lbs)  -  3.38'  •  (12,694  lbs)  -  Mheei  =  0 
Mheei  =  38,694  ft-lbs 

M 

Mu=1.7Hf(D  +  L)  M„=  — ^ 

9 

=  1.7  •  1.3  •  (38,694  ft-lbs)  ^  85,514  ft-lbs 

=  85,514  ft-lbs 

=  95,016  ft-lbs  =  1 140  in-kips 

Because  the  dimensions  (i.e.,  d)  for  both  the  stem  and  heel  and  M„  are 
essentially  the  same,  the  reinforcement  used  for  the  stem  is  also  used  for  the  heel 
(i.e.,  #1 1  @  9"  c-c  =>  As  =  2.08  in^  per  ft  of  wall). 
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Minimum  Effective  Depth: 


j  2.4956  •  M„ 

i  b 

2.4956-  1140  in-kips 
— 

=  15.4"  <  d  okay 

Shear  Capacity  of  Heel: 

V„=1.7-  (D  +  L) 

=  1.7  ■  (18,000  lbs  +  2400  lbs  -  12,694  lbs) 

=  13,100  lbs 


(j)*Vc  ^shear  *2*  >/  f  *c  b*d 

=  0.85-2-  ^/iOOOpsT  12"-19.5" 

=  9';  ISO  ihc 


(|)-V,>1.3-(Vu-(l)-Vc) 

1 .3-(V„  -  (j.- Vc)  =  1 .3-(  1 3, 1 00  lbs  -  25, 1 59  lbs) 
<0  okay 


A.3.4  Shear  capacity  of  the  toe 

The  critical  section  for  shear  capacity  in  the  toe  is  at  a  distance  d  from  the 
interface  of  the  toe  and  the  stem,  where  d  =  19.5".  Because  the  length  of  the  toe 
is  short,  moment  capacity  is  not  checked. 


V,oe=  4293  lbs -4 13  lbs 
=  3880  lbs 

V„=1.7-  (D  +  L) 

=  1.7  -  (3880  lbs) 

=  6596  lbs 


1.375  ft 
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<l)‘Ve  ^shear  '2'  \j  f'c  b‘d 


fV,>1.3-(Vu-(|)-Vc) 


=  0.85  •  2  •  74000  psi  12"-19.5" 
=  25,159  lbs 


1.3-(V„-  fVc)=  1.3-(6596  lbs-25,159  lbs) 
<0  okay 


Use  #1 1  @  9"  c-c  =>  As  =  2.08  in^  per  ft  of  wall 

Figure  A-5a  shows  the  steel  reinforcing  detailing,  determined  previously 
(note:  development  lengths  need  to  be  checked),  wherein  the  stem,  heel,  and  toe 
were  treated  as  singly  reinforced  members  (i.e.,  reinforcement  is  only  on  one  face 
of  the  member).  However,  good  engineering  judgment  would  dictate  that  the 
reinforcing  steel  in  the  toe  and  heel  be  continuous  for  the  length  of  the  base,  such 
as  shown  in  Figure  A-5b.  By  doing  so,  any  potential  issues  regarding  develop¬ 
ment  length  of  the  base  reinforcing  are  avoided,  and  additionally,  the  base  would 
better  withstand  load  reversals  that  might  occur  during  the  life  of  the  wall. 


Figure  A-5.  Proposed  steel  reinforcement  detailing.  A  minimum  4"  cover  is  required 

Notes:  The  primary  purpose  of  the  structural  design  computations  presented 
previously  was  to  ensure  that  the  wall  dimensioned  for  global  stability  was  also 
feasible  from  a  structural  design  standpoint.  The  structural  design  is  not  com¬ 
plete,  and  no  attempts  were  made  to  optimize  the  design  (i.e.,  refine  the  spacing 
of  the  reinforcing  steel,  or  reduce  the  reinforcement  in  the  upper  portions  of  the 
stem). 
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Table  A-1 

Conversion  Factors,  Non-SI  to  SI  Units  of  Measurement  Used  in  This 
Appendix 

Multiply 

By 

To  Obtain 

feet 

0.3048 

meters 

foot-pounds 

1.355818 

joules 

kips  per  square  inch 

6.894757 

megapascals 

Inches 

25.4 

millimeters 

inch-kips 

112.9848 

newton-meters 

pounds  (force)  per  square  foot 

47.88026 

pascals 

pounds  (force)  per  square  inch 

0.006894757 

megapascals 

pounds  (mass) 

0.4535924 

kilograms 

pounds  (mass)  per  cubic  foot 

16.01846 

kilograms  per  cubic  meter 

square  inches 

645.16 

square  millimeters 
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Appendix  B 

Notation,  Sign  Convention,  and 
Earth  Pressure  Expressions 


The  notation  shown  in  Figure  B-1  is  used  throughout  this  report.  All  the 
variables  shown  in  this  figure  are  presented  in  their  positive  orientation.  Addi¬ 
tionally,  expressions  for  the  classical  Mononobe-Okabe  active  and  passive 
dynamic  earth  pressures  are  presented  (e.g.,  Ebeling  and  Morrison  1992,  Chap¬ 
ter  4),'  as  well  as  expressions  for  the  slope  of  the  corresponding  failure  planes. 

Notation: 

Qh  -  horizontal  acceleration 
Gv  =  vertical  acceleration 
g  =  acceleration  due  to  gravity 
H  =  wall  height 
Kae  =  active  thrust  coefficient 
KpE  =  passive  resistance  coefficient 
kh  =  horizontal  inertial  coefficient 
ky  ==  vertical  inertial  coefficient 
Pae  =  active  thrust  of  soil 
PpE  =  passive  resistance  of  soil 
W  =  weight  of  failure  wedge 

*  References  cited  in  this  appendix  are  included  in  the  References  section  at  the  end  of 
the  main  text. 
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Figure  B-1.  Sign  Convention.  All  quantities  are  shown  in  their  positive  orientation  (adapted 
from  Davies,  Richards,  and  Chen  1986) 


Wa  -  weight  of  active  failure  wedge 

fVp  =  weight  of  passive  failure  wedge 

aAE  =  inclination  to  horizontal  of  failure  plane,  active  case 

apE  =  inclination  to  horizontal  of  failure  plane,  passive  case 

Pa  =  slope  of  backfill,  active  case 

Pp  =  slope  of  backfill,  passive  case 

5a  =  interface  friction  angle,  active  case 
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5p  =  interface  friction  angle,  passive  case 
Yt  =  total  unit  weight  of  soil 
Qa  -  wall  batter  angle,  active  case 

0p  =  wall  batter  angle,  passive  case 

(|)  =  angle  of  internal  friction  of  soil 

V|/  =  seismic  inertia  angle 

The  following  expressions  are  the  classical  Mononobe-Okabe  active  and 
passive  dynamic  earth  pressure  expressions  (e.g.,  Ebeling  and  Morrison  1992). 

Mononobe  and  Matsuo  (1929)  and  Okabe  (1926)  derived  these  expressions  by 
modifying  Coulomb's  classical  earth  pressure  expression  to  account  for  seismic- 
ally  induced  horizontal  and  vertical  inertia  forces  in  the  assumed  failure  wedge. 

a.  Active  case: 

(B-1) 


cos\(l>-y/-ej 


cos{y/)-cos^  {6^)-cos{y/  +  0^  +SJ- 


1  + 


I  cos  (,S^+y/  +  0J-  cos  -  0^ ) 


(B-2) 


b.  Passive  case: 

(B-3) 


cos\(j>-y/  +  0;f 


cos {y/) ■  cos" (0  )-cos(ii/-0+S^ 


1- 


I  cos  (S+y^-0J-  cos  (0  -  0^ ) 


(B-4) 


The  expressions  for  the  angles  of  inclination  of  the  failure  planes  with  the 
-  horizontal  originally  derived  by  Zarrabi-Kashani  (1979),  active  case,  and  Ebeling 
and  Morrison  (1992),  passive  case,  and  presented  in  Ebeling  and  Morrison 
(1992)  as  Equations  37  and  61,  respectively,  were  determined  to  give  erroneous 
results  for  <])  <  v|/  +  0.  Equations  37  and  61  in  Ebeling  and  Morrison  (1992)  are 
valid  when  the  wall  batter  angle  is  zero  (i.e.,  0  =  0),  and  for  nonzero  batter  angles 
for  most  levels  of  shaking  (i.e.,  (|)  >  \|/  +  0).  The  following  expressions  for  active 
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and  passive  slip  plane  angles  were  derived  by  Michalowski'  and  numerically 
verified  by  the  authors. 

a.  Active  case: 


^ae  =  tan 


sin  +  y/)  •  c^p  +  cos(^  +  S^+e^+y/) 
cos(y0„  ■t-y/)-c^p-sm{(l>  +  d„+e^+y/) 


(B-5a) 


c 


AE  ~ 


sin(^  +  <5’J-cos(6>^+<5’^+ty) 
sin  - 1^)  •  cos  ((9„  - /?J 


(B-5b) 


b.  Passive  case: 


,  sin (yg^ - y/) •  Cpp  +  cos {(f>  +  Sp-e^  +  y/) 
cos(,pp-y/)-Cpi;+s\n{<l)  +  dp-ep+y/) 


(B-6a) 


Cpp  — 


sin(<Z>  +  )  •  cos(^p  -Sp-xjf) 

sin(^^  +  Pp-y/)-  cosiOp  -  J3^ ) 


(B-6b) 


^  R.  Michalowski,  2002,  Personal  communication  regarding  the  derivation  of  the  angle  of 
inclination  of  active  and  passive  failure  planes  under  seismic  loading.  Department  of  Civil 
Engineering,  University  of  Michigan,  Ann  Arbor,  MI, 
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Appendix  C 

Displacement-Controlled 
Design  Procedure 


The  displacement-controlled  approach  for  design  of  earth  retaining  structures 
was  originally  proposed  by  Richards  and  Elms  (1979)*  and  is  further  developed 
in  Ebeling  and  Morrison  (1992),  Section  6.3.  The  approach  builds  off  the 
Newmark  sliding  block  analysis  procedure  and  provides  a  method  of  assessing 
the  global  stability  of  retaining  walls  against  sliding  and  overturning  based  on  the 
explicit  choice  of  an  allowable  permanent  displacement. 

During  a  seismic  event,  the  maximum  inertial  load  that  can  be  sustained  by 
the  structural  wedge  (in  the  case  of  a  cantilever  retaining  wall)  is  limited  by  the 
maximum  transmissible  acceleration  (N*-g).  N*-g  is  the  value  of  acceleration 
(with  coefficient  N*)  imparted  to  the  driving  wedge  and  structural  wedge  that 
results  in  a  factor  of  safety  against  sliding  equal  to  1 .0.  Recall  that  the  cantilever 
retaining  wall,  contained  within  the  structural  wedge,  was  sized  for  the  normal, 
static  load  case  using  a  factor  of  safety  against  sliding  equal  to  1 .5  (calculations 
given  in  Appendix  A).  So  the  maximum  transmissible  acceleration  is  a  value  of 
acceleration,  N*-g,  when  imparted  to  the  center  of  mass  of  both  the  structural 
wedge  and  the  driving  wedge,  that  reduces  the  factor  of  safety  against  sliding 
from  its  static  value  of  1.5  down  to  a  value  of  1.0.  The  value  of  g  is  a 
function  of  the  geometry  of  the  cantilever  retaining  wall,  the  structural  and 
driving  wedges,  unit  weights  of  the  reinforced  concrete  structural  members,  as 
well  as  the  unit  weight(s)  of  the  soil(s),  and  the  shear  strengths  of  the  backfill,  the 
foundation,  and  the  structure-to-foundation  interface.  Note  that  the  value  of  N*-g 
is  not  a  function  of  the  design  ground  motion.  Also,  consideration  is  given  only 
to  N*-g  corresponding  to  the  movement  of  the  wall  away  from  the  backfill,  as 
this  is  generally  less  than  N*-g  corresponding  to  the  movement  of  the  wall  toward 
the  backfill. 

The  underlying  premise  of  the  displacement-controlled  procedure  is  pre¬ 
sented  in  Figure  C-1.  As  illustrated  in  Figure  C-1  at  time  a,  when  the  horizontal 
inertial  coefficient  kh  =  N*,  sliding  occurs.  Relative  movement  of  the  wall  and 


*  References  cited  in  this  appendix  are  included  in  the  References  section  at  the  end  of 
the  main  text. 
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Figure  C-1.  Incremental  failure  by  base  sliding  (adapted  from  Richards  and  Elms  1979) 

soil  continues  until  their  velocities  become  equal,  as  occurs  at  time  b  in 
Figure  C-l.  The  wall  incrementally  displaces  away  from  the  backfill  each  time 
horizontal  acceleration  Uh  > 

The  Newmark  sliding  block  procedure  is  implemented  in  the  computer 
program  CWROTATE  (Ebeling  and  White,  in  preparation),  which  was  used  to 
evaluate  candidate  ground  motions  for  use  in  the  FLAG  analyses,  the  selection 
criterion  being  that  enough  permanent  relative  displacement  occurs  that  active 
earth  pressure  conditions  exist. 

C.1  CWROTATE  Analysis  of  Cantilever  Retaining 
Wall 

In  order  to  perform  a  Newmark  sliding  block  analysis,  A^-g  is  required. 
CWROTATE  computes  iV**g  by  performing  an  equilibrium  analysis  of  the 
structural  wedge.  The  forces  acting  on  the  structural  wedge  are  shown  in 
Figure  C-2. 

Summing  the  forces  in  the  horizontal  direction  results  in: 

^  ~  ^AE,  heel  “  ~  ^  (C“la) 

where 

r  =  A7'*tan(^z>;)  (C-lb) 
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Structural  Wedge 


fVw  =  weight  of  structural  wedge 

k/,  =  horizontal  inertial  coefficient 

=  vertical  inertial  coefficient 
=  0  {assumed) 

T=  base  shear  reaction  force 
N'  =  base  normal  reaction  force 
5  =  interface  friction  angle 


Figure  C-2.  Forces  acting  on  the  structural  wedge 


N'  =  W^  for  k^=0  (C-lc) 

where  Pae,  heel  is  the  resultant  force  of  the  static  and  dynamic  stresses  acting  on 
the  vertical  section  through  the  heel  and  =  effective  internal  friction  angle  of 

foundation  soil  (use  of  5^,  i.e.,  interface  friction  angle  between  base  and 
foundation  soil,  may  be  also  appropriate). 

Substituting  Equations  C-lb  and  C-lc  into  Equation  C-la  and  solving  for  h 
gives: 


(C-2) 

W 

In  order  to  solve  Equation  C-2  PAEMei  is  needed.  The  formulation  in 
C  WROTATE  uses  the  (active)  wedge  method  of  analysis  and  is  based  on  the 
engineering  formulation  given  in  Appendix  A  of  Ebeling  and  Morrison  (1992)  to 
determine  the  value  for  PAEMet-  A  sweep  search  algorithm  was  devised  and 
implemented  within  CWROTATE  to  ensure  that  the  thrust  force  computed  for  a 
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particular  wedge  geometry  by  the  wedge  method  corresponds  to  the  maximum 
force  (i.e.,  the  pseudo-static  active  earth  pressure  force,  ^AEMet)-  In  this  wedge 
method  of  analysis  a  planar  slip  surface  is  assumed.  The  sweep  search  investi¬ 
gates  the  potential  slip  planes  for  the  trial  wedges  by  CWROTATE  in  1-degree 
increments.  Bilinear  backfill  surfaces  may  be  analyzed  via  CWROTATE  (i.e.,  a 
sloping  backfill  surface  up  to  an  elevation  at  which  the  backfill  surface  is 
horizontal).  For  CWROTATE,  Ebeling  and  Morrison  (1992)  Appendix  A, 
formulation  for  the  (active)  wedge  method  of  analysis  was  expanded  to  consider 
both  effective  stress  (c',  (j)')  and  total  stress  shear  strength  parameters  for  the 
soil.  The  current  2002  version  of  CWROTATE  is  capable  of  analyzing  “diy” 

(i.e.,  moist)  backfills  with  no  water  table  within  the  backfill  nor  foundation. 

Based  on  the  engineering  formulations  given  in  Chapter  4  and  Appendix  A  of 
Ebeling  and  Morrison  (1992),  the  computer  program  CWROTATE  is  being 
expanded  by  Ebeling  to  include  consideration  of  “wet”  sites;  i.e.,  partially 
submerged  backfills  and  submerged  foundations. 

For  the  cantilever  wall  retaining  dry,  level  backfill  that  was  designed  using 
the  procedures  outlined  in  Appendix  A,  CWROTATE  computed  N*g  equal  to 
0.22g.  This  value  of  maximum  transmissible  acceleration  was  also  obtained  using 
Equation  C-2  in  conjunction  with  the  Mononobe-Okabe  equation  for  ’Paemsi- 
Note  that  because  the  Mononobe-Okabe  relationship  (given  in  Appendix  B  of 
this  report  and  in  Chapter  4  of  Ebeling  and  Morrison  1992)  is  a  function  of  kh,  an 
iterative  process  is  used  to  compute  the  kh  {-N*)  value  for  this  wall  by  this 
alternative  method  of  analysis.  The  value  for  PAEMeei  computed  by  CWROTATE 
was  also  verified  using  the  Mononobe-Okabe  relationship.  Recall  that  the 
Mononobe-Okabe  formulation  is  valid  for  homogenous  backfills  with  constant 
backfill  surface  slopes  (including  level  backfills  as  in  this  problem).  Thus,  the 
wedge  method  of  analysis  as  implemented  in  CWROTATE  is  applicable  to  more 
types  of  Corps  retaining  wall  configurations  than  is  the  Mononobe-Okabe 
equation. 

In  order  to  evaluate  the  candidate  motion,  a  site  response  analysis  was 
performed  on  the  profile  shown  in  Figure  3-2  using  a  modified  version  of 
SHAKE91  (Idriss  and  Sun  1992).  The  acceleration  time-history  SG3351  was 
specified  as  an  outcrop  motion  in  a  SHAKE  analysis,  and  the  interlayer  ground 
motion  was  computed  at  a  depth  corresponding  to  approximately  midheight  of 
the  cantilever  retaining  wall.  This  motion  was  used  in  conjunction  with  the  N*g 
in  CWROTATE  to  perform  a  Newmark  sliding  block-type  analysis  of  the 
cantilever  retaining  wall.  The  results  of  the  analysis  are  shown  in  Figure  C-3. 

As  may  be  observed  from  this  figure,  a  permanent  relative  displacement  of 
2.147  in.  (54.5  mm)  is  predicted,  which  is  sufficient  to  ensure  active  earth 
pressure  conditions  (minimum  displacement  for  (seismic)  active  earth  pressures 
behind  the  20-ft-  (6-m-)  high  retaining  wall  was  determined  using  the  Ebeling 
and  Morrison  (1992)  Table  1  data). 
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C.2  Displacement-Controlled  Design  Procedure 

The  following  is  a  summary  of  the  steps  listed  in  Ebeling  and  Morrison 
(1992),  Section  6.3.2,  for  the  displacement-controlled  design  procedure  for  a  wall 
retaining  dry  backfill,  for  which  ky  is  assumed  equal  to  zero.  The  procedure 
pertains  to  the  global  stability  of  the  wall  and  is  repeated  here  for  reference 
purposes  only.  It  is  applicable  to  walls  whose  geometry  has  already  been 
established  (e.g.,  an  existing  wall  or  established  during  a  static  design  step). 
These  steps  are  geared  toward  hand  computations  with  the  earthquake  hazard 
defined  in  terms  of  peak  acceleration  and  peak  velocity  values  only  and  not  for 
time-history  analysis  of  permanent  wall  displacement  (as  is  the  case  when  using 
CWROTATE  to  analyze  retaining  walls).  Interface  friction  between  the  driving 
wedge  and  the  structural  wedge  is  assumed  to  be  zero  for  illustrative  purposes 
only  in  the  following  equations.  The  reader  is  directed  to  Ebeling  and  Morrison 
(1992)  for  more  detailed  information  on  the  simplified  engineering  procedure  and 
for  the  case  when  interface  friction  between  the  structural  wedge  and  the  driving 
wedge  is  nonzero. 

a.  Determine  the  value  for  the  average  site-specific  peak  horizontal 
acceleration  A  -g  and  the  value  for  the  average  peak  horizontal  velocity  V 
at  the  site. 

b.  Calculate  the  maximum  transmissible  acceleration  N*-g  coefficient  N*. 
An  iterative  method  consisting  of  the  following  steps  is  required  to 
compute  the  value  for 

(1)  Using  an  assumed  value  for  A^,  compute  the  value  for  the  dynamic 
active  earth  pressure  force  PAEMei  using  the  Mononobe-Okabe 
relationship  given  in  Appendix  B  (or  in  Chapter  4  of  Ebeling  and 
Morrison  1992). 

(2)  Calculate  the  values  of  the  shear  force  required  for  equilibrium  along 
the  base  of  the  wall,  T: 

T^PAEMet-^^-^h  (C-3) 

(3)  Calculate  the  normal  force,  N: 

N  =  W  (C-4) 

(4)  Calculate  the  value  for  the  ultimate  shear  force  along  the  base  of  the 
wall  Tull  using  the  smaller  of  the  following  two  values  for  Tun: 

T^„=N-taniS,)  (C-5) 

or 

r„„=)V-tan((zi',)  (C-6) 
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where 


5b  =  interface  friction  angle  between  base  of  the  wall 
and  foundation  soil 

=  effective  internal  friction  angle  of  foundation  soil 

(5)  If  the  value  for  T  is  not  equal  to  the  value  for  Tuu,  adjust  the  value 
used  for  N*  and  repeat  steps  (1)  through  (4)  until  T  =  The 
resulting  value  for  N*  is  equal  to  the  meiximum  transmissible 
acceleration  coefficient. 

c.  Calculate  the  permanent  relative  displacement  dr  using  Whitman  and 
Liao  (1985): 

,  r495-F' 

where 

A-g  =  expressed  in  units  of  inVsec^ 

V  =  expressed  in  units  of  in./sec 
dr  =  expressed  in  units  of  inches 
g  =  386  in/sec^ 

d.  The  value  for  dr  must  be  consistent  with  those  movements  that  are 
required  to  develop  the  dynamic  active  earth  pressure  used  in  Step  b(l). 
Refer  to  Table  1  in  Section  2.2.2  of  Ebeling  and  Morrison  (1992). 

e.  Dimension  the  wall  so  that  the  Corps’  overturning  criterion 
(Headquarters,  U.S.  Army  Corps  of  Engineers,  1989)  is  satisfied. 


exp 
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Appendix  D 

Specifying  Ground  Motions  in 
FLAG 


As  briefly  outlined  in  Chapter  3,  dynamic  analyses  can  be  performed  with 
FLAG,  wherein  user-specified  acceleration,  velocity,  stress,  or  force  time- 
histories  can  be  input  as  an  exterior  boundary  condition  or  as  an  interior 
excitation.  A  parametric  study  was  performed  to  determine  the  best  way  to 
specify  the  ground  motions  in  FLAG  for  earthquake  analyses.  The  parametric 
study  involved  performing  a  series  of  one-dimensional  site  response  analyses 
using  consistently  generated  acceleration,  velocity,  and  stress  time-histories. 
Generally,  earthquake  ground  motions  are  not  defined  in  terms  of  force  time- 
histories  and  therefore  were  not  considered  in  the  parametric  study. 

Using  the  acceleration  time-history  SG3351,  a  SHAKE  analysis  was 
performed  on  a  5  percent  damped,  nondegrading  profile,  shown  in  Figure  D-1, 
wherein  the  SG3351  was  specified  as  an  outcrop  motion.  Interlayer  acceleration 
and  stress  time-histories  were  computed  at  bedrock  and  at  depths  of  50, 35,  and 
25  ft  (15, 10,  and  8  m),  in  addition  to  the  surface  acceleration  time-history  being 
computed.  Interlayer  velocity  time-histories  were  computed  by  integrating  the 
interlayer  acceleration  time-histories  using  the  trapezoidal  rule.  The  interlayer 
acceleration,  velocity,  and  stress  time-histories  were  used  as  base  motions  in  a 
series  of  FLAG  analyses,  in  which  the  acceleration  time-histories  at  the  surface  of 
the  FLAG  profiles  were  computed.  The  profiles  used  in  the  FLAG  analyses  were 
comparable  to  the  225-ft  (69-m)  SHAKE  profile  down  to  the  depths 
corresponding  to  the  interlayer  motions.  This  is  illustrated  in  Figure  D-1 .  An 
elastic  constitutive  relation,  with  5  percent  Rayleigh  damping,  was  used  to  model 
the  soil  layers  in  the  FLAG  profiles.  The  central  frequency  of  the  damping  was 
set  to  the  fundamental  frequency  of  the  respective  FLAG  profiles. 

Fourier  amplitude  spectra  (FAS)  and  5  percent  damped,  pseudo-acceleration 
spectra  (PSA)  were  computed  from  the  acceleration  time-histories  of  the  surface 
motions  of  the  SHAKE  and  FLAG  profiles.  Error  analyses  were  performed  on 
the  spectra  corresponding  to  the  different  profiles  and  different  types  of  specified 
input  motions.  In  the  error  analyses,  the  spectra  for  the  SHAKE  motions  were 
used  as  the  “correct”  motions.  The  word  “correct”  does  not  imply  that  SHAKE 
correctly  models  the  behavior  of  an  actual  soil  profile  subjected  to  earthquake 
motions.  Rather,  SHAKE  gives  the  analytically  correct  motion  for  a  visco-elastic 
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profile  with  constant  damping  applied  to  all  frequencies  of  motion.  On  the  other 
hand,  the  FLAG  models  used  in  this  study  give  only  numerical  approximations  of 
the  correct  analytical  solution.  The  percent  error  in  the  FLAG  spectral  values,  as 
a  function  of  frequency,  was  computed  using  Equation  D-1 : 


„  FLAG  Spectral  Value -SHAKE  Spectral  Value 

Error  (%)  = - - - - - x  1 00% 

SHAKE  Spectral  Value 


(D-1) 


In  Figures  D-2a,  D-3a,  D-4a,  and  D-5a,  the  5  percent  damped  PSA  are  shown 
for  the  surface  motions  of  the  respective  profiles  and  types  of  input  motions.  The 
errors  in  the  PSA  computed  using  Equation  D-1  are  presented  in  Figures  D-2b, 
D-3b,  D-4b,  and  D-5b.  Because  the  scales  of  these  figures  do  not  allow  a  clear 
presentation  of  the  errors  in  the  spectral  values  for  the  motions  corresponding  to 
acceleration  and  velocity  input  motions.  Figures  D-2c,  D-3c,  D-4c,  and  D-5c  are 
provided.  In  these  figures,  only  the  errors  in  the  spectral  values  for  the  motions 
corresponding  to  acceleration  and  velocity  input  motions  are  plotted. 

As  may  be  observed  from  these  figures,  the  specification  of  the  input  motion 
in  terms  of  stress  time-histories  gives  the  least  accurate  results.  Although  it 
appears  that  specification  of  input  motions  in  terms  of  acceleration  time-histories 
gives  slightly  more  accurate  results  than  using  velocity  time-histories,  the 
apparent  error  may  be  related  to  how  the  interlayer  velocity  time-histories  were 
computed  from  the  interlayer  acceleration  time-histories.  Regardless,  the  errors 
corresponding  to  specifying  the  motions  in  terms  of  acceleration  and  velocity 
time-histories  are  approximately  equal.  The  input  motions  used  in  the  FLAG 
analyses  of  the  retaining  walls  were  specified  in  terms  of  acceleration  time- 
histories. 

Similar  error  analyses  were  performed  on  the  FAS  of  the  surface  motions, 
with  similar  conclusions  being  drawn.  However,  the  explicit  solution  scheme 
used  in  FLAG  introduces  considerable  high-frequency  noise,  which  dominated 
the  scales  of  the  error  plots  of  the  FAS.  Accordingly,  these  plots  are  not 
presented. 


D2 


Appendix  D  Specifying  Ground  Motions  in  FLAG 


Shear  Wave 

Layer 

SHAKE  Profile  FLAG  Profiles 

Velocity 

(fps) 

Thickness 

(ft) 

225ft  225ft  50ft  35ft 

25  ft 

_.J25 _ 

_ A.  . 

"  ■iiuiiiiLU  —xm 

600 

6 

_ 

6S0 

4 

6 

- 

-.im. _ 

.  .5.  .  . 

855 

10 

''• - interlayer  motion - 

BASE 

930 

15 

V - interlayer  motion - ^  BASE 

1015 

20 

interlayer  motion  BASE 

1100 

25 

1185 

30 

1250 

30 

1325 

35 

1385 

35 

HALFSPACE  \  interlayer  f  EASE 
motion 

Figure  D-1 .  Profiles  used  in  the  parametric  study  to  determine  the  best  way  to  specify  earthquake 
ground  motions  in  FLAG 
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Figure  D-2.  Results  of  the  error  analysis  for  the  225-ft  (69-m)  profile:  (a)  Pseudo-acceleration 
spectra  of  surface  motions,  (b)  percent  error  in  pseudo-spectral  values 
corresponding  to  input  motions  specified  in  terms  of  stress,  acceleration,  and 
velocity  time-histories,  (c)  percent  error  in  pseudo-spectral  values  corresponding 
to  input  motions  specified  in  terms  of  acceleration  and  velocity  time-histories  (i.e., 
same  as  (b)  except  plotted  using  a  different  scale) 
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Figure  D-3.  Results  of  the  error  analysis  for  the  50-ft  (15-m)  profile;  (a)  Pseudo-acceleration 
spectra  of  surface  motions,  (b)  percent  error  in  pseudo-spectral  values 
corresponding  to  input  motions  specified  in  terms  of  stress,  acceleration,  and 
velocity  time-histories,  (c)  percent  error  in  pseudo-spectral  values  corresponding  to 
input  motions  specified  in  terms  of  acceleration  and  velocity  time-histories  (i.e., 
same  as  (b)  except  plotted  using  a  different  scale) 
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Period  (sec) 


Figure  D-4.  Results  of  the  error  analysis  for  the  35-ft  (1 1-m)  profile:  (a)  Pseudo-acceleration 
spectra  of  surface  motions,  (b)  percent  error  in  pseudo-spectral  values 
corresponding  to  input  motions  specified  in  terms  of  stress,  acceleration,  and 
velocity  time-histories,  (c)  percent  error  in  pseudo-spectral  values  corresponding 
to  input  motions  specified  in  terms  of  acceleration  and  velocity  time-histories  (i.e., 
same  as  (b)  except  plotted  using  a  different  scale) 
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Error  ( percent)  Error  (  percent)  Pseudo-Spectral  Acceleration  (g) 


Figure  D-5.  Results  of  the  error  analysis  for  the  25-ft  (8-m)  profile:  (a)  Pseudo-acceleration 
spectra  of  surface  motions,  (b)  percent  error  in  pseudo-spectral  values 
corresponding  to  input  motions  specified  in  terms  of  stress,  acceleration,  and 
velocity  time-histories,  (c)  percent  error  in  pseudo-spectral  values  corresponding 
to  input  motions  specified  in  terms  of  acceleration  and  velocity  time-histories  (i.e. 
same  as  (b)  except  plotted  using  a  different  scale) 
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ah  Horizontal  acceleration 

Ag  Cross-sectional  area 

Ck  Interface  stiffness  ratio 

dr  Relative  displacement 

Ec  Elastic  modulus 

/  Frequency 

fmax  Highest  frequency 

/J  Compressive  strength  of  concrete 

Fsr  Lateral  seismic  force  component 

Fv  Vertical  component  of  the  resultant  force  acting  on  the  heel  section 

Fh  Horizontal  component  of  the  resultant  force  acting  on  the  heel 
section 

Fly  Force  acting  on  the  top  node  of  element  /  and  at  time  increment  j 

pbottom  acting  on  the  bottom  node  of  element  i  and  at  time 

increment  j 

FDM  Finite  difference  method 

FEM  Finite  element  method 

FLAG  Fast  Lagrangian  Analysis  of  Continua 
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g  Acceleration  due  to  gravity 

G  Shear  modulus 

h/  Length  of  beam  element  / 

H  Height  of  the  wall 

/  Second  moment  of  area,  or  moment  of  inertia 

kh  Dimensionless  horizontal  inertial  coefficient 

kn  Normal  stiffness  for  interface  element 

ks  Shear  stiffness  for  interface  element 

ky,j  Vertical  inertial  coefficient  at  time  increment  j  (assumed  to  be  zero) 

K  Lateral  earth  pressure  coefficient 

K'  Bulk  modulus 

Ka  Active  earth  pressure  coefficient 

Kae  Active  dynamic  earth  pressure  coefficient 

Ki  Dimensionless  interface  stiffness  number  for  initial  loading 

KpE  Passive  dynamic  earth  pressure  coefficient 

Ko  At-rest  earth  pressure  coefficient 

Ksi  Initial  shear  stiffness  of  the  interface 

Kurj  Unload-reload  stiffness  number  for  interfaces 

m  Seismically  induced  bending  moment  on  a  section  of  the  stem 
portion  of  a  cantilever  wall 

M  Magnitude  of  the  earthquake 

M  Mean  magnitude  of  the  earthquake 

MDE  Maximum  Design  Earthquake 

Hj  Dimensionless  stiffness  exponent 

N'  Normal  reaction  force  to  forces  acting  on  the  structural  wedge 
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N*-g  Maximum  transmissible  acceleration 
OBE  Operational  Basis  Earthquake 
pga  Peak  ground  acceleration 

Pa  Atmospheric  pressure 

Paemci  Resultant  force  of  the  static  and  dynamic  stresses  acting  on  the 
vertical  section  through  the  heel 

PE.stem  Resultant  force  of  the  static  and  dynamic  stresses  acting  on  the  stem 

of  the  wall 

Pheei  Resultant  forces  acting  on  the  heel  section 

Pj  Total  force  acting  on  the  stem  or  heel  section  at  time  increment  j 

Pstatic  Total  resultant  force  prior  to  shaking 

Psiem  Resultant  forces  acting  on  the  stem  section 

PSHA  Probabilistic  seismic  hazard  analyses 

R  Site-to-source  distance  of  the  ground  motion 

R  Mean  site-to-source  distance  of  the  ground  motion 

Rjj  Failure  ratio  for  interface 

s  Seismically  induced  shear  moment  on  a  section  of  the  stem  portion 

of  a  cantilever  wall 

S  Slider;  yield  strength  of  interface  element 

SASW  Spectral  analysis  of  surface  waves 

T  Base  shear  reaction  force  to  the  static  and  dynamic  stresses  acting 
on  the  vertical  section  through  the  heel 

T'  Tensile  strength  of  interface  element 

Shear  wave  velocity 

WUS  Western  United  States 

y,  Vertical  distance  from  the  base  of  the  retaining  wall  to  the  center  of 

beam  element  i 
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Y  Vertical  distance  from  the  base  of  the  retaining  wall  to  the  resultant 
force 

Yheei  Vertical  distance  from  the  base  of  the  wall  to  the  point  of 
application  of  the  resultant  forces  acting  on  the  heel  section 

Yj  Vertical  distance  from  the  base  of  the  retaining  wall  to  the  point  of 
application  of  the  total  resultant  force  acting  on  the  stem  or  heel 
section  at  time  increment  j 

Ysiaiic  Vertical  distance  from  the  base  of  the  retaining  wall  to  the  point  of 
application  of  the  total  resultant  force  acting  on  the  stem  or  heel 
section  prior  to  the  shaking  (i.e.,  Yj  at  j  =  0) 

Ysiem  Vertical  distance  from  the  base  of  the  wall  to  the  point  of 

application  of  the  resultant  forces  acting  on  the  stem  section 

(YP)hcei  Resultant  forces  multiplied  by  the  corresponding  vertical  distance 
above  the  base  at  which  they  act  for  the  heel  section 

(Y-P)siem  Resultant  forces  multiplied  by  the  corresponding  vertical  distances 
above  the  base  at  which  they  act  for  the  stem  section 

a  Mass-proportional  damping  constant 

P  Stiffness-proportional  damping  constant 

y  Unit  weight  of  the  soil 

y,  Total  unit  weight  of  the  soil 

jw  Unit  weight  of  water 

8  Interface  friction  angle 

A/  Length  of  an  element 

APj  Incremental  dynamic  forces  at  time  increment  j 

AY  Vertical  distance  from  the  base  of  the  retaining  wall  to  the  point  at 
which  APj  acts 

Azm/„  Smallest  width  of  a  zone  in  the  normal  direction  of  the  interfacing 
surface 

X  Wavelength 

^  Critical  damping  ratio 
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V 


Poisson's  ratio 


p  Mass  density 

(7 1  j  Average  lateral  stress  acting  on  element  i  and  at  time  increment^ 

On  Normal  stress  acting  on  the  interface,  and  determined  iteratively  in 
FLAG  by  first  assuming  a  small  value  for  4  and  then  constructing 
the  wall 

^bottom  gtress  acting  on  the  bottom  of  element  i  and  at  time 

increment  j 

<jfj  Lateral  stress  acting  on  the  top  of  element  i  and  at  time  incrementy 

(j)  Internal  friction  angle 

CO  Angular  frequency  associated  with  ^ 
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